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ABSTRACT 

In this work nonperturbative aspects of quantum gravity are investigated using the lattice 
formulation, and some new results are presented for critical exponents, amplitudes and invariant 
correlation functions. Values for the universal scaling dimensions are compared with other nonper¬ 
turbative approaches to gravity in four dimensions, and specifically to the conjectured value for the 
universal critical exponent = 1/3. It is found that the lattice results are generally consistent with 
gravitational anti-screening, which would imply a slow increase in the strength of the gravitational 
coupling with distance, and here detailed estimates for exponents and amplitudes characterizing 
this slow rise are presented. Furthermore, it is shown that in the lattice approach (as for gauge 
theories) the quantum theory is highly constrained, and eventually by virtue of scaling depends 
on a rather small set of physical parameters. Arguments are given in support of the statement 
that the fundamental reference scale for the growth of the gravitational coupling G with distance 
is represented by the observed scaled cosmological constant A, which in gravity acts as an effective 
nonperturbative infrared cutoff. In this nonperturbative vacuum condensate picture a fundamen¬ 
tal relationship emerges between the scale characterizing the running of G at large distances, the 
macroscopic scale for the curvature as described by the observed cosmological constant, and the 
behavior of invariant gravitational correlation functions at large distances. Overall, the lattice 
results suggest that the slow infrared growth of G with distance should become observable only 
on very large distance scales, comparable to A. It is hoped that future high precision satellite 
experiments will possibly come within reach of this small quantum correction, as suggested by a 
vacuum condensate picture of quantum gravity. 
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1 Introduction 


In this work nonperturbative aspects of the ground state for quantum gravity will be discussed, 
based on the lattice theory. So far the lattice formulation represents the only known first-principle 
method for reproducing correctly the low energy properties of non-Abelian gauge theories, including 
confinement and chiral symmetry breaking. It is hoped therefore that a lattice approach to quantum 
gravity will shed useful light on the low energy properties of quantum gravity as well. A key 
aspect of this method is a recognition of the importance of Wilson’s modern interpretation of 
the renormalization group mm as it applies to perturbatively non-renormalizable theories mm, 
including gravity. In previous work the elegant lattice formulation of gravity of Regge and Wheeler 
mm was used to compute a number of observable quantities expected to be relevant for the ground 
state properties of quantum gravity (for a detailed discussion of the Feynman path integral approach 
to quantum gravity the reader is referred to 0)- The lattice formulation implies the existence of 
an ultraviolet cutoff, nevertheless in the real world such a cutoff would presumably arise from short 
distance details derived from an underlying more fundamental theory, such as higher derivative 
gravity, supergravity or string theory. Nevertheless, it is expected that a softer cutoff would lead 
to significant short distance modifications of gravity, while leaving the quantum infrared behavior 
largely unchanged. It is these universal long distance effects that form the subject of the present 
paper. As in QED and QCD, it will turn out that these quantum infrared modifications to gravity 
are intrinsically non-local. 

In four dimensions (and for the Euclidean theory) it was found that for gravity two phases are 
possible, a pathological gravitational screening phase for G < Gc, and an anti-screening phase for 
G > Gc- It has been known for some time that the screening phase corresponds to a branched 
polymer, with no physically acceptable continuum limit, and the lattice results were therefore 
interpreted as suggesting that ultimately the only physical acceptable phase is the strong coupling 
phase for G > Gc- Eurthermore, it was found that in this phase the average local curvature 
approaches zero towards the critical point at Gc, indicating that in this phase the recovery of the 
semiclassical limit for gravity appears to be possible. In previous work detailed estimates were 
given for the location of the critical point at Gc, for the critical exponents and scaling dimensions 
characterizing the growth of the gravitational coupling with distance, and for the scaling behavior 
of gravitational correlation functions. Of central importance in these results is the value for the 
universal critical exponent v, related to the derivative of the beta function for G at the hxed 
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point. Here more refined estimates for the critical point and scaling dimensions will be provided; 
the analysis will later be extended to correlation functions of invariant operators at fixed geodesic 
distances. In this context the present discussion includes both local operators, as well as extended 
ones such as the Wilson loop and the correlation between gravitational Wilson loops. In previous 
work it was argued that the gravitational Wilson loop provides information about the macroscopic 
curvature, and therefore about the hoped-for recovery of the semiclassical limit. This last result is 
quite different from what is found in gauge theories, since in gravity the gravitational Wilson loop 
has no significance for the static potential. Also, it was shown earlier that the the area law for the 
gravitational Wilson loop provides a connection between the nonperturbative scale ^ that arises in 
nonperturbative gravity and the macroscopic large-scale curvature, and thus the observed effective 
cosmological constant A. Here it will be argued that the numerical results so far are consistent 
with many of the previous answers, including the conjecture that the exponent v is exactly equal 
to one third in four dimensions. The latter part of the paper will therefore deal with a detailed 
discussion of the possible physical significance of having an exponent v exactly equal to one third in 
four dimensions, as this relates to a number of physical consequences, such as the scale dependence 
of G and the behavior of invariant gravitational correlation functions at large separation. 

The structure of the paper is as follows. In Section 2 the form of the discretized lattice gravi¬ 
tational Feynman path integral will be recalled, and basic notation will be established. Here the 
basic dehnitions for local averages and their fluctuations will be laid out as well. Section 3 will 
introduce basic diffeomorphism invariant correlation functions, and show how these can be tran¬ 
scribed to the lattice theory. Section 4 will extend the previous discussion to correlation functions 
involving operators that are not necessarily local, such as the correlation between smeared oper¬ 
ators, the dehnition of the gravitational Wilson loop, and the form and basic expected properties 
of correlations between these loop operators. Section 5 will recall the basic fundamental scaling 
assumptions for the gravitational path integral, and how those assumptions and definitions affect 
the critical behavior of various averages, fluctuations and correlations dehned in the previous sec¬ 
tions. Section 6 summarize how the quantum continuum limit in lattice quantum gravity should 
be taken, in accordance with the general principles of the renormalization group. A discussion is 
provided to show how the interplay between the bare coupling constants, the critical point and the 
correlation length ^ lead to a definite expression for the running of Newton’s G. In addition, it is 
shown how the prediction for a running of G can be described in terms of universal quantities, to 
leading order in the vicinity of that fixed point, and specihcally in terms of universal exponents 
and amplitudes. Sections 7,8,9 and 10 later provide details on the how the numerical calculations 
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are performed, and on the methods by which the universal critical exponents and amplitudes are 
extracted from the numerical results. A discussion is given to show the overall consistency of the 
results, based on a variety of different observables and methods of analysis. At the end of Section 
10 the results obtained from a variety of different observables are compared in a comprehensive 
summary table. Two additional tables later provide a comparison between the lattice results for 
the universal exponents and the values obtained by other nonperturbative methods. In Section 10 
a separate comparison table is provided for four dimensions, and a second table is added later for 
the case of three dimensions. It is argued that the numerical results so far are consistent with the 
expectation that the universal critical exponent v for quantum gravity in four dimensions is equal 
to one third. Section 11 then discusses the physical implications of having an exponent v exactly 
equal to one third as it applies to various local averages, fluctuations and correlation functions 
introduced earlier in the paper. It is shown that many lattice results become particularly simple 
and perhaps more transparent for this choice of exponent. The numerical calculation also supply 
a value for several critical amplitudes which appear in the running of G in the vicinity of the fixed 
point at Gc- Section 12 at the end of the paper is devoted to a discussion of the curvature corre¬ 
lation function, and how this correlation can, in suitable cases via the field equations, be related 
to the analogous correlation function for matter density fluctuations. The final section presents 
some conclusions, and elaborates further on a suggestive analogy between the vacuum condensate 
picture derived from lattice quantum gravity and the well understood nonperturbative properties 
of non-Abelian gauge theories, and specifically the case of lattice QGD. 


2 Path Integral, Invariant Local Gravitational Averages and their 
Fluctuations 


In this section the basic definitions for diffeomorphism invariant gravitational averages and corre¬ 
lations will be recalled briefly, in a form suitable for later discussions. Here the starting point for 
a nonperturbative formulation of quantum gravity is the discretized form for the Feynman path 
integral for pure gravity [8], written originally as 

Zc= dn[g^u] ex.p{-I[g^^]} , ( 1 ) 
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with invariant gravitational action 

= I (ao - I ^ + • • •) (2) 

and DeWitt invariant functional measure [9] 

f Mapu] = /n (Vfi'(^)) n ■ (3) 

In the above expression k~^ = SttG with G the bare Newton’s constant, Aq the bare cosmological 
constant, and qq a possible higher derivative coupling m- In the absence of matter fields the 
DeWitt invariant measure for pure gravity in four dimensions corresponds to the simple choice 
a = 0. In the following we will only consider the case oq = 0, i.e. no higher derivative i?^-type 
terms. 

The continuum Feynman path integral given above is generally ill-defined, and has to be for¬ 
mulated more precisely by introducing a suitable discretization m- The last step is particularly 
important for nonperturbative calculations, where the nontrivial invariant measure over the 
plays a key role. Regge and Wheeler proposed an elegant discretization of the classical gravitatio¬ 
nal action EE], which forms the basis for the lattice formulation of quantum gravity discussed 
in this paper. Once the measure and the path integral have been discretized, the ultimate goal 
then becomes to recover the original continuum theory of Eq. © in the limit of a small lattice 
spacing (this limit is rather subtle, and involves in a nontrivial way fundamental aspects of the 
renormalization group). This approach then leads, as a starting point, to the following discrete 
form for the Euclidean Feynman path integral for pure gravity 

= j d/i[;2]exp{-/[/2]} , (5) 

with lattice gravitational action 

= 6hAh + a blAl/Vh + • • • ) ( 6 ) 

h 

^ A well known problem of the Euclidean path integral formulation is the conformal instability of the classical 
gravitational action [111 112]. The latter is seen by considering conformal transformations with Q a 

positive function. Then the Einstein-Hilbert action transforms into 

I[5l = {D.'^R + 6g'"''d^ndM) ■ (4) 

which can be made arbitrarily negative by choosing a rapidly varying conformal factor fl. The wrong sign for the 
kinetic term of the Q fields then implies that the Euclidean gravitational functional integral is possibly badly divergent, 
depending on the detailed nature of the gravitational measure contribution dfi[g^^] (the “entropy” or phase space 
part), and specifically its behavior in the regime of strong fields and rapidly varying conformal factors. 
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and lattice functional measure 


/ /*oo 

d^^f] = / n n diij 




( 7 ) 


Here the sum over hinges h in four dimensions corresponds to a sum over all lattice triangles with 
area A^, with deficit angles describing the curvature around them. 

In the discrete formulation a functional integration over metric is replaced by an integration 
over squared edge lengths, which are taken as fundamental variables in the discrete theory. The 
basis for this step is a rather direct correspondence between the squared edge lengths in a four- 
simplex and the induced metric within that same simplex. Within each n-simplex s one can define 
a metric in terms of unit vectors e* pointing along the edges 


9 ij{s) = ei-ej , ( 8 ) 

with 1 < i,j < n, and a positive definite quantity in the Euclidean case. In terms of the edge 
lengths lij = I e* — Cj I, or conversely 


gijis) — ^ {^oi + ^oj ~ ^ij) (9) 

for a simplex based at 0. This last result then provides a key connection between the metric gp,y{x) 
in the continuum and the lattice degrees of freedom if, which is essential in establishing a fairly 
unambiguous relationship between lattice and continuum operators, just as is the case in ordinary 
lattice gauge theories. It is also known that the lattice action in Eq. ([5] generally reduces to the 
continuum one of Eq. ([T]for smooth enough field configurations [IT], and that it contains the correct 
physical degrees of freedom for gravity in the weak field limit, namely transverse-traceless (massless 
spin two) modes [T8] . 

The general aim of the calculations presented later will be to evaluate the lattice path integral 
exactly by numerical means, by performing a (correctly weighted) sum over all fields configurations, 
without relying on the weak field expansion, or an expansion around suitable saddle points or some 
other approximate scheme, which generally tends to involve a number of assumptions on what 
configurations (smooth or otherwise) might or might not play a dominant role in the path integral 
(indeed the general expectation for such path integrals is that smooth field configurations tend to 

® In the following we will deal almost exclnsively, as is customary in lattice field theories, with dimensionless 
quantities. Thus the couplings Ao and G appearing in the continuum theory will be expressed from the start in 
units of the fundamental lattice ultraviolet cutoff A — 1/a [14] . As is standard procedure in ordinary lattice field 
theories and lattice gauge theories [151116|. the latter is then later set equal to one, which means that from then 
on all observable quantities, correlations and couplings are measured in units of this fundamental ultraviolet cutoff. 
The actual value for the ultraviolet cutoff (in MeV or cm~^) is later determined by comparing suitable physical 
quantities, see Eqs. (Ton and Ton towards the end of the paper. 
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have measure zero). Here the functional integration over edge lengths is highly nontrivial, due to 
the constraint coming from the generalized triangle inequalities [expressed in the function Q[lfj] 
in Eq. ([7])], which is placed there in the Euclidean formulation to insure that all edge lengths, 
triangle areas, tetrahedra and simplex volumes are strictly positive. The discrete gravitational 
measure in Zl of Eq. ([5|) can then be regarded as a regularized version of the DeWitt continuum 
functional measure [9]. Also, a bare cosmological constant term is essential for the convergence of 
the path integral, while curvature squared terms allow one to further control the fluctuations in 
the curvature [HI IS]. It is generally understood that these last terms are generated by radiative 
corrections within a perturbative diagrammatic treatment in the continuum. In practice, and for 
obvious phenomenological reasons, one is nevertheless only interested eventually in a limit where 
effective higher derivative contributions are negligible compared to the rest of the action, oq — 0. 

In this limit the theory depends, in the absence of matter and after a suitable rescaling of the 
metric (in the continuum) or the edge lengths (on the lattice), only on one bare parameter, the 
dimensionless coupling k/y/\Q. Indeed already in the continuum one finds in d dimensions under a 
rescaling of the metric 

9iiu = UJ , ( 10 ) 

with 00 a constant, that the cosmological constant term Xo^/g turns into Aq so that a 

subsequent rescaling 

G ^ Ao ^ Ao (11) 

leaves only the dimensionless combination unchanged. Clearly only the latter combination 

has physical meaning in pure gravity, and in particular one can always choose the scale cu = Ag 
so as to adjust the volume term to have a unit coefficient. Equivalently, this shows that it seems 
physically meaningless to discuss separately the renormalization properties of G and Aq. Without 
any loss of generality one can therefore set the bare cosmological constant Aq = 1 in units of 
the ultraviolet cutoff [TT|. The latter contribution controls the overall scale for the edge lengths, 
contains (like a mass term) no derivatives in the continuum, and does not affect the construction 
of a suitable lattice continuum limit, which is determined by the relative interplay between the 
curvature and volume terms. It seems therefore redundant to vary Aq, as this will only change 
the overall length scale, without any discernible effect on the quantum lattice continuum limit. In 
the continuum a similar result can be derived; there on can show that the renormalization of Aq is 
gauge- and scheme-dependent, and that only the renormalization of G is independent of the choice 
of gauge condition [22l [23l IH] . 
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Some partial information about the behavior of physical correlations can be obtained indi¬ 
rectly from averages of suitable local invariant operators. In [19] a set of diffeomorphism invariant 
gravitational observables, such as the average curvature and its fluctuation, were introduced. Ap¬ 
propriate lattice analogs of these quantities are easily written down, making use of the following 
well understood correspondences 

Vdix) ^ Vh 

hinges hDx 

^R{x) Y 

hinges hD x 

V9 ^ 4 ^ {5hAhf/Vh ■ ( 12 ) 

hinges hDx 

An overall numerical normalization coefficient has been omitted on the r.h.s., since it will depend 
on how many hinges are actually included in the summation. In the following we will not consider 
any further higher derivative terms, which means that the subsequent discussion will be limited 
almost exclusively to the first and second type of operators. 

First consider the average local curvature, defined as 

< f d'^x y/g R{x) > 

< f d^Xy/g> 

The above quantity is relevant for parallel transports of vectors around an elementary, infinitesimal 
parallel transport loop, and is by construction manifestly diffeomorphism invariant. On the lattice 
one first notes that it is preferable to define quantities in such a way that variations in the average 
lattice spacing /q ~ y/<R> are compensated by a suitable multiplicative factor, determined entirely 
from dimensional considerations. It would be possible to adjust Aq in Eq. ([5|) to achieve Iq = 1, 
but here we choose to have simply Aq = 1 in units of the ultraviolet cutoff A = I/a. As stated 
previously, in the following all quantities will be expressed in units of this fundamental cutoff a, 
whose value, as is customary in lattice gauge theories, is set initially equal to one, a = 1. In the 
case of the average local curvature a useful lattice definition is therefore [laiiii 

W(t) . (14) 

<2^h^h> 

Note that by construction this quantity is dimensionless, and consequently if all edge lengths are 
rescaled by a common factor it remains unchanged. Again, this choice factors out an entirely 
irrelevant overall length scale (a phenomenon peculiar to gravity, which does not arise in ordinary 
lattice gauge theories). 




8 





A second quantity of interest is the local curvature fluctuation 


Xn{k) ~ 


< (f d^Xy/gRY > — < f d'^x y/g R 
< f d^Xy/g> 


A suitable lattice transcription of this last quantity is 


( 15 ) 


Xn{k) 


< 

<EhVh> 


(16) 


Note that in the functional integral formulation of Eqs. o and ([5]) both the average curvature Tl{k) 
and its fluctuation XTl{k) can be obtained by taking derivatives of the functional Zl in Eq. ([5|) 
with respect to k. Therefore on the lattice one has 


n{k) ~ 


1 


d 


<V> dk 


InZL 


(17) 


and 

X7e(A:) ~ ^^^InZi , (18) 

just as the analogous continuum quantities in Eqs. (fTH)) and (fTCD can be obtained as derivatives of 
the expression in Eq. ([T]). In a similar way, the average volume per site is defined as 


<'">-• («) 

h 

and again one has 

Furthermore, its fluctuation xv can also be obtained as a second derivative of Z^ with respect to 
the bare cosmological constant Aq. A simple scaling argument, based on neglecting the effects of 
curvature terms entirely (which vanish in the vicinity of the critical point), is found to give a rather 
accurate estimate for the average volume per edge 


< V; > ~ 


2(1 + ad) 
Ao d 


1 

—y - 

d=4, a=0 2 Aq 


( 21 ) 


In four dimensions numerical simulations agree quite well with this simple formula. Finally, a set 
of exact sume rules can be derived from the scaling properties of the action and measure in Eq. ([5]). 
As an example, for the case of the dl"^ measure one finds the following exact lattice Ward identity 


2Ao < ^ Eh > - k<Y^ ShAh > - Ni = 0 , (22) 

h h 

which is easily derived from Eq. ([S]) and the definitions in Eqs. (|17l) and (I20p . Here A^o represents 
the number of sites in the lattice, and the averages are defined per site. For the hypercubic lattices 
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used in this paper, A^i = 15 A^O; ^2 = 50 A'q, = 36 Nq and iV 4 = 24 A^q- The above exact identity 
can be a useful tool in establishing the numerical convergence of the integration method used for 
the lattice path integral. @ 


3 Diffeomorphism Invariant Gravitational Correlation Functions 


Generally in a quantum theory of gravity the physical distance between any two points x and y in 
a fixed background geometry is determined by the metric 

d{x,y\g) = min f . (23) 

S J t{x) 

Because of quantum fluctuations the latter depends on the metric or, equivalently, in the lattice 
case on the edge length configuration considered. Correlation functions of local operators need 
to account for this fluctuating distance, and as a result these correlations have to be computed 
at some fixed geodesic distance between a set of given spacetime points mills]. In addition, in 
gravity one generally requires that the local operators involved should be coordinate scalars. In 
principle one could also smear such operators over a small region of spacetime, an option which 
will be discussed later. It is also possible to compute nonlocal gravitational observables in analogy 
to what is done in Yang-Mills theories, by defining objects such as the gravitational Wilson loop 
(which carries information about the parallel transport of vectors around large loops, and therefore 
about large scale curvature) [26l ET] EH] ES], or the correlation between Wilson lines closed by 
the lattice periodicity (which can be used for extracting the static potential in quantum gravity) 
isniEi]. One more different type of gravitational correlation was studied in |32j . 

A fundamental correlation function in te quantum theory of gravity is the one associated with 
the scalar curvature, with physical points x and y separated by a fixed geodesic distance d 

Gnid) < y/g R{x) y/g R{y) 5{\x - y\-d) >c ■ (24) 

It is then straightforward to define the same type of object on the lattice. If the lattice deficit 
angles are averaged over a number of contiguous hinges which share a common vertex, one is lead 
to consider the connected lattice correlation function at fixed geodesic distance d 

Gnid) = <'^2 6hAh 2 6h'Ah'6{\x - y\ - d) >c . (25) 

hZ)x h'Dy 

^ As an example, in practice one can achieve that the l.h.s. of Eq. (I22t is zero to about one part in 10® for 200fc 
individual lattice edge length configurations containing around 25 million simplices. 
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The need to compute physical distances between points for any given metric (or edge length) 
field configuration complicates the problem considerably, as compared for example to ordinary 
gauge theories, where the distance between points is assigned a priori based on a fixed immutable 
underlying lattice structure. 

For the curvature correlation at fixed geodesic distance one expects at short distances (i.e. 
distances much shorter than the gravitational correlation length a power law decay 

< ^ R{x) y/g R{y) 6{\x-y\-d) >c ~ , (26) 

with the power characterized by a universal exponent n; how n is related to another calculable 
universal critical exponent (in particular u) will be discussed further below. 

One notes on the other hand that for sufficiently strong coupling (large G, or small k) fluctu¬ 
ations in different spacetime regions largely decouple (the kinetic or derivative term in Eqs. ([T]) or 
([^ is responsible for coupling fluctuations in different regions, and it comes with a coefficient 1/G). 
In this regime one then expects a faster, exponential decay, controlled by the correlation length ^ 


<y/gR{x)^R{y)5i\x-y\-d)>c e . (27) 

d 3> 5 


This last result shows that the fundamental gravitational correlation length if nonzero, can be 
dehned thro^h the long-distance decay of the connected invariant correlations at fixed geodesic 
distance d. Ij Note also that the behavior in Eq. ([26]) is expected to hold at short distances, 
i.e. distances much larger than the fundamental lattice spacing but significantly shorter than the 
correlation length, Iq d whereas the behavior in Eq. (|27p is expected to hold at much larger 
distances, d S> ^ ^ /q- In either case, in order to reach the lattice continuum limit the distances 


considered need to be much larger than the fundamental lattice spacing, d,^ 3> Iq- This last 

^Practical useful methods for calculating such diffeomorphism invariant correlations were described in detail in 
|25 |. For each given metric configuration which, properly weighted, contributes to the path integral one needs to 
compute both the geodesic distance between any two points, as well as the correlation between a set of given invariant 
operators centered at those points. By far the most time consuming part of the calculation is the determination of 
the actual physical distance between any two given points for an assigned background metric configuration. The 
latter part can be done by generating a large number of random walks that start at one of the two points, and 
then obtaining the physical distance from the shortest walk. Alternatively, the geodesic distance can be determined 
directly from the propagator, and specifically the exponential decay in distance of a covariantly coupled lattice scalar 
field propagator with a given mass. Either way, the calculation is then later repeated for every metric configuration 
contributing to the chosen ensemble, resulting eventually in the sought-after final average. 

® This rather general result can be proven easily by using the same type of arguments used in ordinary field 
theories (and lattice gauge theories) to show that Euclidean correlation functions generally decay exponentially at 
strong coupling. There one shows that it takes n actions of the kinetic (or hopping) term to connect, via the shortest 
possible lattice path, two points that are n lattice sites apart. A similar result holds for lattice gravity, where the 
relevant kinetic or hopping term is the curvature {R) contribution, proportional to 1/G [33]. Of course in the extreme 
limit of infinite G, due to the absence of a kinetic term, fluctuations in the fields at different spacetime locations 
completely decouple, and in this limit the correlation length shrinks to zero (or more precisely, to one lattice spacing). 
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constraint is referred to as the scaling limit, where short distance lattice artifacts are presumably 
washed out, and the true (and physically relevant) continuum limit is expected to emerge. Later 
it will be shown, from rather elementary scaling considerations, that the exponent n in Eq. (|26p is 
related to the so-called correlation length exponent u in four dimensions by n = 4 — l/v. 

Another key result of relevance here lies in the fact that the local curvature fluctuation defined 
in Eqs. (flSl) and (fT6]l is directly related to the connected curvature correlation of Eqs. ([241) and 
(ESI) at zero momentum 


f d^xf d^y < ^g(x)R(x) g{y) R{y) >e 
< f d^xy/g{xj > 


(28) 


a well-known and rather useful result already in ordinary field theories. This connection will be 
used extensively further below, and its relevance will lie in the fact that it allows one to relate the 
exponent v, obtained for example from the curvature fluctuation of Eq. m, to the physical corre¬ 
lation in Eqs. (I24|) and ()26l) . The latter can in turn be related to other physical correlations, such as 
the matter density correlation, by the use for example of the effective, long distance gravitational 
field equations. 


4 Correlations Between Smeared and Nonlocal Operators 


The discussion up to this point has dealt with local operators and their correlations, i.e. operators 
defined at a point x in spacetime, or on the lattice at a single lattice point. Some rather mild 
nonlocality does in fact appear, due to the circumstance that the gravitational action involves, 
via the affine connection and the Riemann tensor, the parallel transport of a test vector around 
an infinitesimally small loop. The latter is encoded on the lattice by the deficit angles, which 
describe the parallel transport of a vector around a loop whose size is comparable to the lattice 
spacing, or in physical terms of size comparable to the ultraviolet cutoff or the Planck length. It is 
nevertheless possible to define smeared operators, which involve a new length scale: the linear size 
of the smearing volume r*. On the basis of rather general renormalization group arguments one 
expects correlations for these operators to have milder short distance divergences. 

Consider first the average of an operators over a spherically shaped smearing region Q(x,rs), 
centered at the point x and of linear size r^. In other words, all points within a physical distance 
Xg from the point in question are considered; on a lattice of course the number of points within a 
given physical neighborhood of the point x with linear size Xg will in general be finite. Then define 
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Figure 1: Correlation function for two smeared operators of size separated by a geodesic 
distance d. 


the smeared operator Os{x) by the spacetime average 

Os{x) = [ d^z^Oiz) . (29) 

J n(a;,rs) 

A natural candidate operator for smearing is of course the scalar curvature, but various curvature 
squared terms would also be viable, for example. A suitable invariant correlation function is then 
defined as 

G,, {d) = < Osix) Os{y) 5i\x -y\-d)>,, (30) 


where again the correlation between the two operators Osix) is taken at a fixed geodesic distance 
d. The general expectation is that the short distance (d > r^) behavior for this correlation function 
is less singular than for correlations of operators dehned at a single point. Nevertheless at larger 
distances d 3> the asymptotic decay of the correlation function should be the same as for the 
one in Eq. (1261) . provided the operators in question have the same quantum numbers. 

A second class of invariant correlation functions for smeared operators involves the correlation 
of parallel transport loops |28t I26j. First note that infinitesimal transport loops appear already in 
the dehnition of the correlation function for the local scalar curvature, as in Eqs. (|24ll and (|25p . 
Next consider the parallel transport of a vector around a loop C which is not infinitesimal; in the 
following this loop will be assumed to be close to planar, a well defined geometric construction 
described in detail in [29]. First dehne the total rotation matrix U(C') along the path C via a 
path-ordered (V) exponential of the integral of the affine connection E^j^ 


U^uiC) 


V exp 


[jjxdx^ 



(31) 
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The lattice action itself already contains contributions from infinitesimal loops, but more generally 
one might want to consider near-planar, but noninfinitesimal, lattice closed loops C. Along such a 
closed loop the overall rotation matrix is given by a product of elementary rotations defined along 
the lattice path 


u^uic) = [n 

s cc 


U. 


s,s+l 


(32) 


In analogy with the inhnitesimal loop case, one expects for the overall rotation matrix 


U^uiC) 




- M 

- U 


(33) 


where oj^y{C) is an area bivector perpendicular to the loop and 5{C) the corresponding deficit 
angle. This will work if the loop is close to planar, so that can be taken to be approximately 
constant along the path C, or defined by some suitable average over the loop. Here by a near-planar 
loop around the point P what is meant is a loop that is constructed by drawing outgoing geodesics 
on a plane through P, so that this unit bivector plays the role of a normal to the loop. A coordinate 
scalar can be defined by contracting the above rotation matrix U(C) with the appropriate unit 
length bivector, namely 

Wc = uj^C) P^^(C) (34) 


where the bivector uJai 3 {C) is taken to be representative of the overall geometric features of the 
loop. Now if the parallel transport loop in question is centered at the point x, then one can define 
the operator Wc{x) by 

Wc{x) = a;^,(C,x) [/^"(C,®) (35) 


with the near-planar loop centered at x and of linear size rc- A suitable invariant correlation 
two-point function for these operators is then defined as 


Gc{d) = < Wc{x) Wc{y) (5(|x -y\-d) >c , (36) 

where again the correlation between the loop operators Wc{x) is taken at some given fixed geodesic 
distance d. Of course for infinitesimal loops one recovers the expressions given earlier in Eqs. (I24p 
and (l2^ . 

In general one needs to specify the relative orientation of the two loops. So, for example, one can 
take the first loop in a plane perpendicular to the direction associated with the geodesic connecting 
the two points, and the same for the second loop; the parallel transport of a vector along this 
geodesic will then be sufficient to establish the relative orientation of the two loops. Nevertheless if 
one is interested in the analog (for large loops) of the scalar curvature, then it will be adequate to 


14 





w 



w 


Figure 2: Correlation function of two infinitesimal parallel transport loops, separated by a geodesic 
distance d. 


perform a weighted sum over all possible loop orientations at both ends. This is in fact precisely 
what is done for infinitesimal loops of size rc ~ a, if one looks carefully at the way the Regge 
lattice action is originally defined. Again here the expectation is that the short distance d > rc 
behavior of this correlation function for extended loop objects is less singular than for correlations 
of operators defined at a point; nevertheless at larger distances d such that rc ^ d ^ the decay 


of these correlation function should be the same as the local ones in Eq. (I26p . 

It is possible to give a more quantitative description for the behavior of the loop-loop correlation 
function given in Eq. (j36D . at least in the strong coupling limit. The following estimate is based 
on the previous results and definitions, and the important analogy and correspondence of lattice 
gravity to non-Abelian gauge theories outlined in |33l I29j . Eirst it will be assumed here that the 
two (near planar) loops are of comparable shape and size, with overall linear sizes rc ~ T and 
perimeter P ~ 2ttL. In addition, the two loops will be separated by a distance d ^ L, and for 
both loops it will be assumed that this separation is much larger than the lattice spacing, d ^ a 
and L ^ a. Then to get a nonvanishing correlation in the strong coupling, large G limit it will be 
necessary to completely tile a tube connecting the two loops, due to the area law arising from the 
use of the Haar measure for the local rotation matrices at strong coupling, again as discussed in 
detail in [29]. In this last paper extensive use is made of a modified hrst order formalism for the 
Regge lattice theory, based on the work of [34], which then allows the separation of metric degrees 
of freedom into local Lorentz rotations and tetrads, as is done in the continuum. Consequently in 
this limit one obtains an area law 



(37) 


Consistency of the above expression with the result for small (infinitesimal) loops given in Eqs. (1261) 
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Figure 3: Correlation function for two large parallel transport loops of size and orientation 
separated by a geodesic distance d. 

and ([271) . and the area law for large loops requires the following limits for the quantity ^o{L) 

Co{L) ~ a and ^o{L) ~ ^ . (38) 

L a L a 

From these results one concludes that the asymptotic decay of correlations for large loops is fun¬ 
damentally different in form as compared to the decay of correlations for infinitesimal loops, with 
an additional factor of ^ appearing for large loops. In other words, the results of Eqs. (1261) and 
(j27p only apply to infinitesimal loops which probe the parallel transport on infinitesimal (cutoff) 
scales, and these results will have to be suitably amended when much larger loops, of semiclassical 
significance, are considered. 


5 Renormalization Group Scaling Relations for Gravity 


In this section some of the basic scaling relations for quantum gravity will be summarized. It is 
by now established wisdom, at least in most field theories besides gravity, that standard scaling 
arguments allow one to determine the scaling behavior of local averages, correlation functions 
and even suitable nonlocal observables such as the Wilson loop from the knowledge of the basic 
renormalization group behavior, and specifically from the universal critical exponents. The latter 
generally characterize the singular behavior of local averages in the vicinity of the critical point. 
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a nontrivial fixed point of the renormalization group (RG) in field theory language. For extensive 
reviews on the subject see for example [EIESKIIIIMIET]. There is by now a rather well established 
body of knowledge in quantum field theory and statistical field theory on this subject, and there is 
no apparent reason why its basic tenets should not apply to gravity as well, with quantum gravity 
describing the unique theory of a massless spin two particle coupled to a covariantly conserved 
energy momentum tensor [8]. 

It is also understood that in the vicinity of a critical point (seen as equivalent to a nontrivial 
fixed point of the renormalization group) long range correlations arise due to the appearance of a 
massless particle. In statistical field theory language, the presence of a massless particle is reflected 
in a divergent correlation length ^ = 1/m, or equivalently a power law in the relevant correlation 
functions. Let us summarize here the basis for the scaling assumptions for local averages, fluctu¬ 
ations and their correlations. ^ In brief, since logZ in either Eq. © or © is both dimensionless 
and extensive, for a volume R ~ it has to have the form 

log Z(k) = fa (0" + fs (I)" (39) 

where here fa and fs are nonsingular functions of dimensionless parameters, is the fundamental 
correlation length (the distance over which quantum fluctuations are strongly correlated), and 
/q ~ a the fundamental lattice spacing or ultraviolet cutoff. Then the free energy or generating 
function, dehned as 

F{k) = -y log Z(k) , (40) 

is expected, based on purely dimensional grounds, to aquire a singular part Fsing{k) such that [Hj 

Fsing{k) ~ r'" • (41) 


If one sets for the nonperturbative correlation length ^ 

m A^ika-k)-\ (42) 

k^kc 

^ It is well established that for theories with a nontrivial ultraviolet fixed point mm, the long distance (and 
thus infrared) universal scaling properties are uniquely determined, up to subleading corrections to exponents and 
scaling amplitudes, by the (generally nontrivial) scaling dimensions obtained via renormalization group methods 
in the vicinity of an ultraviolet fixed point [m m 133 EHi EZ]- These sets of results form the basis of universal 
predictions for, as an example, the (perturbatively nonrenormalizable) nonlinear sigma model |38l I39| . The latter 
provides today one of the most accurate test of quantum field theory |40| . after the g — 2 prediction for QED (for a 
comprehensive set of references, see [HE], and the references therein). It is also a well established fact of modern 
renormalization group theory that in lattice QCD the scaling behavior of the theory in the vicinity of the asymptotic 
freedom ultraviolet fixed point unambiguously determines the universal nonperturbative scaling properties of the 
theory m, as quantified by physical observables such as hadron masses, vacuum condensates, decay amplitudes, the 
QCD string tension etc. |421I43 |. 
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where is the correlation length amplitude, kc the critical point and u the correlation length 
exponent characterizing the divergence of ^ at the critical point, then one obtains for the singular 
part of the free energy 

Fsing ~ {kc-kf-. (43) 

k — ^kc 

One concludes that a divergent correlation length signals the presence of a phase transition, and 
this in turn leads to the appearance of nonanaliticities in thermodynamic quantities such as Z{k) 
and the free energy F{k). The origin of these nonanaliticities in Z{k) lies therefore in the divergence 
of ^ in the vicinity of the critical point at kc, where the theory becomes scale invariant. 

The following results then follow more or less immediately from the dehnitions in Eqs. © or 
and in Eq. (HZD. Near the singularity the average curvature behaves as 

n{k) -An{kc-kf , (44) 

k—ykc 

with a curvature exponent 6 related to the exponent introduced earlier in Eq. ()42l) by the scaling 
relation 6 = di^ — 1. Consequently the presence of a phase transition can already be inferred 
directly from the appearance of nonanalytic terms in invariant local averages, such as the average 
curvature. 1^ Similarly, one has for the curvature fluctuation defined in Eq. (|16p . using Eqs. (IlSp 
and ()i2]) . 

Xiz(k) ~ 6 Ati (kc - . (45) 

k—>-kc 

Again, scaling [Eqs. (1421) and (1431) ] relates the exponent 6 appearing in the curvature fluctuation 
to ly, so that the exponent in Eq. (I45p is simply 1 — 6 = 2 — dv. These results show that from 
suitable averages one can extract the correlation length exponent v [defined in Eq. (1421) ]. without 
even a need to compute an invariant two-point function, such as the ones in Eqs. (IMP . (IMP . (IMP 
and (|27p . provided scaling holds. Eurthermore, in the vicinity of the critical point kc one can trade 
the distance from the critical point for the correlation length and obtain the following equivalent 
result relating the quantum expectation value of the local curvature to the physical correlation 
length 

no ■ (46) 

k — ykc. 

This last expression is obtained from Eqs. (|42p and (I44p . using 5 = du — 1. Matching of di¬ 
mensionalities here can always be achieved by supplying appropriate powers of the lattice spacing 
Iq ^ a or, equivalently, the Planck length Ip = y/G. 

® An additive constant could be present in Eq. (|44l). but the evidence so far points to this constant being consistent 
with zero for the Regge lattice gravity theory. 
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In addition, the above results allow one to relate the fundamental scaling exponent v of Eq. (1421) 
to the scaling behavior of some correlation functions at large distances. Thus, for example, the 
curvature fluctuation of Eq. (jlSp is related to the connected scalar curvature correlation of Eq. (1261) 
evaluated at zero momentum 




/ / dS < y/gR{x) y/gRjy) >c 

< f > 


r\j 

k — ykc 


Ay- {kc 




(47) 


It follows that a divergence in the curvature fluctuation is indicative of long range correlations, 
corresponding to the presence of a massless particle, the graviton. Close to the critical point one 
expects, in the scaling limit, i.e. for physical distances much larger than the fundamental lattice 
spacing, a power law decay in the geodesic distance d, as in Eq. (I26l) . 


< y/gR{x)y/gR{v) >c 


1 


(48) 


\x-y\^cxy |x-y|2"- 

After inserting the correlation function from Eq. ()48p in the expression of Eq. m and then inte¬ 
grating over a region of size ^ one obtains immediately for the power in Eq 


n = d — l/v . (49) 

Thus knowledge of the scaling exponent ly uniquely determines the power in Eqs. (|48D and (1261) . 

So, the scaling theory for gravity just outlined implies a universal relationship between various 
quantities and the exponents that appear in them. Of course the nonperturbative scaling exponents 
can be determined, in principle, separately for each individual observable. Nevertheless scaling 
theory, based on the assumption of the existence of a massless particle in the vicinity of an ultraviolet 
hxed point, immediately implies a direct (and testable) relationship between the scaling behavior of 
various quantities, such as the ones in Eqs. (fT3]) . (fT^ . ([26]) . (I27P and (1421) . Ultimately, the physical 
relevance of the above results is that Eq. ()42p gives, when solved for k or G, the running of G 
with scale in the vicinity of the fixed point at Gc- Thus, again, Eas. ([41P and (1451) are useful for an 
accurate determination of the universal scaling exponent v, and this quantity in turn determines 
the scaling behavior of invariant curvature correlation functions in Eqs. (1261) and (j48p as a function 
of geodesic distance. 

® Note that in weak field perturbation theory [26] < y/gR{3:)y/gR{y) >c ~ < d^h{x)d^h{y) > ~ l/ja; — y\‘^'^^, 
which is quite different from the result in Eq. (I48II unless u = 2/(d — 2), which is only correct for d close to two, where 
Einstein gravity becomes perturbatively renormalizable and the corrections to free field behavior become small. 
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6 Continuum Limit of Lattice Quantum Gravity 


The long distance behavior of quantum field theories is, to a great extent, determined by the scaling 
behavior of the relevant coupling constants under a change in momentum scale. Asymptotically 
free theories such as QCD lead to vanishing gauge couplings at short distances, while the opposite 
is true for QED. In general the fixed point(s) of the renormalization group need not be at zero 
coupling, but can be located at some finite Gc, leading to nontrivial fixed points or more complex 
limit cycles [aiiiiiiKiB]. 

These general ideas are realized concretely in the analytic 2 + e expansion for gravity, and 
reappear later in essentially the same form in lattice gravity in four dimensions. In the 2 + e 
perturbative expansion for gravity [201 ED ED EH] one analytically continues in the spacetime 
dimension using dimensional regularization, and applies perturbation theory about d = 2, where 
Newton’s constant becomes dimensionless. A similar method is quite successful in determining the 
critical properties of the 0(n)-symmetric nonlinear sigma model above two dimensions [Uj- In the 
framework of this expansion the dimensionful bare coupling is written as Gq = G, where A is 
an ultraviolet cutoff (corresponding on the lattice to a momentum cutoff comparable to the inverse 
average lattice spacing, A~ I/Zq ~ 1/®)- There were originally some known technical difficulties 
with this expansion due to the presence of kinematic singularities for the graviton propagator in two 
dimension (the Einstein action is a topological invariant in d = 2), but these have been overcome 
recently. In addition, one can show that a gauge-choice dependent renormalization of the bare 
cosmological constant Aq can be completely reabsorbed into an overall rescaling of the metric, with 
no physical consequences. A double expansion in G and e = d — 2 then leads in lowest order to a 
gauge-independent nontrivial fixed point in G above two dimensions 

r)G 

with /3o > 0 for pure gravity. To lowest order the ultraviolet fixed point is then at Gc = l//^o(f^“2). 
Integrating Eq. ()50l) close to the nontrivial fixed point one obtains for G > Gc 

m„ = Aexp(-y , (51) 

where mo an integration constant, with dimensions of a mass or inverse length, expected to be 
associated with some physical scale. It is rather natural here to identify this scale with the inverse of 
the gravitational correlation length (^ = m“^), or some equivalent scale associated with the physical 
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large-scale curvature [281129] . Note that the derivative of the beta function at the fixed point defines 
the critical exponent u, which to this order is independent of /3o, l3'{Gc) = —{d — 2) = —l/v. 

The previous results clearly illustrate how the lattice continuum limit should be taken. It 
corresponds to A ^ oo, G ^ Gc with the physical scale ^ = 1/m held constant; thus for fixed 
lattice cutoff the continuum limit is approached by tuning G to Gc- In four dimensions the universal 
critical exponent is defined by [see Eq. ([l2]) ] 

r\G) = m{G) ~ A™A|G(A)-G,r , (52) 

G^Gc 

where A = 1/a is the inverse lattice spacing, and the nonperturbative mass scale m = 1/^ is defined 
as the inverse of the correlation length, with a nonperturbative but calculable amplitude. The 
cutoff independence of the nonperturbative mass scale m implies 

A^m(A,G(A)) = 0 . (53) 

Comparing results in Eqs. ([5T]l and (l5^ one obtains 

/3'(Ge) = -1/iy , (54) 


so that the universal exponent ly is directly related to the derivative of the Callan-Symanzik /3 
function for G in the vicinity of the ultraviolet fixed point. Thus computing v is equivalent to 
computing the universal derivative of the beta function at Gc- 

It is easy to see here that the value of ly determines the running of the effective coupling G(/u) 
in the vicinity of the fixed point, where // is an arbitrary momentum scale. The renormalization 
group tells us that in general the effective coupling will grow or decrease with length scale r = l//i, 
depending on whether G > Gc or G < Gc, respectively. This result follows from the fact that the 
genuinely nonperturbative physical mass parameter m = is itself scale independent, and obeys 
therefore the rather simple Callan-Symanzik renormalization group equation 


/U ^ ^(a^)) = A ^ {Amp\ G(/i) - GcT } = 0 . (55) 

Here again, by virtue of Eq. (|52p . the second expression on the right-hand-side is appropriate in 

the vicinity of the ultraviolet fixed point at Gc- Nevertheless the above discussion is not necessarily 

As a concrete example, in the 2 + e expansion for pure gravity using the background field method one finds at 
two loop order = (d — 2) -|- |(d — 2)^ + 0((d — 2)®) |23 |. Consistency of this expansion generally requires a 
smooth background with a small Aq in Eq. ©• Nevertheless a renormalization of Aq there is later undone by the 
metric rescaling of Eq. (nni. so that the only physical (and gauge-choice independent) running is in the gravitational 
coupling G. 
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limited to just a region in the immediate vicinity of Gc, more generally, if one defines the function 
F{G) via 

= m = AF{G{A)) , (56) 


then, from the usual definition of the Callan-Symanzik /3 function /3(G) = dG{A)/dlogA, one 
obtains 


/3(G) = - 


F{G) 


(57) 


F'{G) ’ 

which shows that the renormalization group /3-function, and thus the running of G(//) with scale, can 
be defined also some distance away from the nontrivial ultraviolet fixed point [including therefore 
the higher order corrections in Eq. (|55ll ]. So, more generally, the running of G(/i) is obtained by 
solving the differential equation 

dn(n\ 

(58) 


dG{ij) 

/X —-- = /3(G(^)) 


dfi 

with /3(G) obtained from Eq. (|57p . hd It is then clear from the previous discussion that the 
physical mass scale m = determines the magnitude of the scaling corrections, and plays a role 
similar to the scaling violation parameter Aj^ in QCD (as in gauge theories, this nonperturbative 
mass scale emerges in spite of the fact that the fundamental gauge boson remains strictly massless 
to all orders in perturbation theory, and consequently does not violate local gauge invariance). 
Eurthermore, as in gauge theories, one expects in gravity that the magnitude of ^ cannot be 
determined perturbatively, and to pin down its value requires a fully nonperturbative approach 
such as the lattice formulation. 

Solving explicitly Eq. ()55p for G{q^), with q an arbitrary wavevector scale, one obtains 


G{q^) = Ge 






(59) 


Here the amplitude of the quantum correction cq is directly related to the constant Am in Eq. (I52p 
by 


Co 


G 


(60) 


One important point is that the magnitude of the quantum correction in Eq. (|59p depends crucially 
on the magnitude of the nonperturbative physical scale Also, the expression in Eq. (l5^ does 
not satisfy general covariance; this in turn can be fixed by performing the replacement —□ 

As an example, for /3(G) = ~^{G — Gc) — b(G — Gc)^, where b is some numerical constant, one obtains for 
the correlation length = m = a Am (G — GcF (l — bv^{G — Gc)^ + 0{{G — Gc)^)), which relates a sub-leading 
correction to /3(G) to the sub-leading correction in m{G). 
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where 'Ci[g^y) is the covariant D’Alembertian for a given background metric g^y{x) 
then leads, from Eq. ([59]), to 


12 


G(n) = Gc 


1 + Co 


I 

0 ^ 


-^2 02 

A set of manifestly covariant effective field equations with a G(n) takes the simple form [46 


This 


(61) 


R^iv - \g^ivR + >^g^iv = Svr G(n) 


(62) 


with the nonlocal contribution coming from the quantum correction in the G(n) of Eq. (|61ll . [)£ 
These nonlocal effective field equations can then be solved for a number of physically relevant 
metrics. For the specific case of a static isotropic metric it is possible to obtain an exact expression 
for G{r) in the limit r ^ 2MG [l6]. The result, for z/ = 1/3 exactly, reads 


G G(r) = G (l + —m^r^ In + ... ) (63) 

\ 3 vr ) 

with m = 1/^. This last result is vaguely reminiscent of the Uehling (vacuum polarization) correc¬ 
tion to the static potential found in QED. Generally the expressions in Eqs. ([59|) . (I6T]) and (l63]) are 
consistent with a gradual slow increase in G with large distance r, and with a modified Newtonian 
potential in the same limit. 

The remainder of this paper will deal therefore with establishing firm values for the nonper- 
turbative amplitudes and exponents defined in the previous sections, and later determining both 
qualitatively and quantitatively their effects on the running of Newton’s G and on the long distance 
behavior of physical correlation functions, such as the ones dehned in the previous sections. 


7 Average Local Curvature 


Next we come to a discussion of the numerical methods employed in this work and the analysis of 
the results. For the reader who is not interested in such details, a separate section later summarizes 

In the lattice theory of gravity only the smooth phase with G > Gc exists (in the sense that an instability 
develops and spacetime collapses onto itself for G < Gc)- This then implies that the gravitational coupling can only 
increase with distance [19]. In other words, a gravitational screening phase does not exist in the lattice theory of 
quantum gravity. This situation appears to be true both for the Euclidean theory in four dimensions and in the 
Lorentzian version in 3 + 1 dimensions |47| . 

If general covariance is to be maintained, then it is virtually impossible here to have a running cosmological 
term in the field equations with a AfD), by virtue of the simple fact that covariant derivatives of the metric vanish 
identically, Va giii> = 0 [24| . 

One can show that this exact solution only exists provided v = l/{d — 1) for d > 4; otherwise no consistent 
solution to the effective nonlocal field equations with GfD) can be found |46| . 
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the most important results obtained so far. As in previous work, the edge lengths are updated by a 
Monte Carlo algorithm, generating eventually an ensemble of configurations distributed according 
to the action and measure of Eq. Details of the method as it applies to pure gravity are 
discussed in OSH], and will not be repeated here. 

In this work lattices of size with L = 4 (256 sites, 3,840 edges and 6,144 simplices), L = 8 
(4,096 sites, 61,440 edges and 98,304 simplices), L = 16 (65,536 sites, 983,040 edges and 1,572,864 
simplices), L = 32 (1,048,576 sites, 15,728,640 edges and 25,165,824 simplices), and L = 64 
(16,777,216 sites, 251,658,240 edges and 402,653,184 simplices) have been considered. These lattices 
are all constructed by conveniently dividing up hypercubes into simplice by introducing suitable 
diagonals m, and periodic boundary conditions are used throughout. In general for a lattice with 
L'^ sites one has 15 edges per vertex and 24 four-simplices per vertex. For these lattices one should 
keep in mind that due to the simplicial nature of the lattice there are many edges per hypercube 
with many interaction terms; as a consequence the statistical fluctuations already for one single 
hypercube can be comparatively small, unless one is very close to a critical point. The results pre¬ 
sented here are still preliminary, and in the future it should be possible to repeat such calculations 
with improved accuracy on even larger lattices. Also, while the overall statistics on the 32^ lattice 
seems adequate, the overall statistics on the 64^^ lattices was yet far too low to be usable in the 
present analysis. 

On the 32^ lattice up to 200,000 consecutive configurations were generated for each value of 
k, and 9 different values for the parameter k were chosen. On the 16^ lattice up to 800,000 
consecutive configurations were generated for each value of k, and 32 different values for k were 
chosen. In addition, results for different values of k can be considered as completely statistically 
uncorrelated, since they originated from unrelated edge length configurations. On the smaller 8^ 
lattice 200,000 consecutive configurations were generated for each value of k. On the 4^ lattice 
two million consecutive configurations were generated for each value of k. To accumulate enough 
statistics, runs were performed around the clock on a 1200 core machine over a period of roughly 
four months. As a result, the increase in accuracy is signihcant compared to the results presented 
in previous work done at the time on a dedicated 32-node cluster [45] . 

In this work the topology is restricted to a four-torus (periodic boundary conditions). One could 
perform similar calculations with other lattices employing different boundary conditions or topology, 
but one expects that the universal long distance scaling properties of the theory to be determined by 
short-distance renormalization effects, which are generally independent of the boundary conditions 
at infinity. A clear example of this is of course the Feynman diagrammatic expansion for gravity in 
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2 + e dimensions, where boundary conditions play no role in the renormalization of the couplings. 
In addition, it will be necessary to impose, based on physical considerations, the constraint that 
the correlation length in lattice units be much larger than the average lattice spacing, and at the 
same time much smaller than the overall linear size of the system, Iq < ^ < Lq, where Lq rsj yi/4 
here is the linear size of the system, and Iq ^ a the (average) lattice spacing. 

As stated earlier, the bare cosmological constant Aq appearing in the gravitational action of 
Eq. ([5]) is set 1 since its value just sets the overall length scale in the problem. The higher deriva¬ 
tive coupling oo was also set to 0 (pure Regge-Einstein action). It is possible to introduce R^-type 
terms in the action, nevertheless in this work these terms were not included in order not to “con¬ 
taminate” the results with the effects of such higher derivative terms. These terms were studied 
extensively in [38], and their effects is generally to stabilize the theory at the expense of nonunitary 
contributions, which cause a visible oscillatory behavior in curvature correlations at short distances. 
The downside of not including any lattice higher derivative terms is that the theory eventually de¬ 
velops instabilities very close to the critical point in G, which need to be handled properly by 
an extrapolation or analytic continuation in G. Nevertheless, as has been shown in [38] and in 
the discussion further below, such an extrapolation or analytic continuation is fairly unambiguous, 
given a large enough amount of high precision numerical data. Indeed such an instability is in 
fact expected on the basis of the well-known Euclidean conformal mode contribution, arising from 
a kinetic energy contribution for the conformal mode with the wrong sign mm Its appearance 
should therefore be regarded as consistent with the full recovery of a continuum behavior in the 
vicinity of the ultraviolet fixed point at Gc- 

For the measure in Eq. ([S]) the above choice of parameters then leads to a well behaved ground 
state for fe < /cc ~ 0.052 for o = 0 [481149] . Given this choice of parameters the system then resides 
in the ‘smooth’ phase, with a fractal dimension close to four; on the other hand for k > kc the local 
curvature can become rather large (‘rough’ phase), and lattice spacetime collapses into a degenerate 
configuration with very long, elongated simplices and thus more akin to a two-dimensional lattice 

[iidHiinKiH]. 

The results obtained for the average curvature TZ [defined in Eq. (fT31) ] as a function of the bare 
coupling k are shown in Figures 4 to 8, on lattices of increasing size with 4^, 8^, 16^ and 32^^ sites. 
Figures 5 and 7 show the 32^^ data by itself. The errors there are quite small, of the order of a 
tenth of a percent or less, and are therefore not visible in the graph. In [35] it was found that as 
k is varied, the average local curvature is negative for sufficiently small k (’smooth’ phase), and 
appears to go to zero continuously at some finite value kc- For k > kc the curvature becomes very 
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large, and the simplices tend to collapse into degenerate configurations with very small volumes 
(< y > / < /^ >^ ~ 0). This collapsed phase corresponds to the region of the usual weak held 
expansion (G ~ 0), characterized by unbounded huctuations in the conformal mode. 

Accurate and reproducible curvature data can only be obtained for k below the instability point 
ku since, as already pointed out in [58], for k > k^ ^ 0.052 an instability develops, presumably 
associated with the unbounded conformal mode. Its signature is typical of a sharp hrst order transi¬ 
tion, beyond which the system tunnels into the rough, elongated phase which is two-dimensional in 
nature with no physically acceptable continuum limit. This instability is caused by the appearance 
of one or more localized singular conhguration, with a spike-like curvature singularity, and is clearly 
driven by the Euclidean Einstein term in the action, and in particular its unbounded conformal 
mode contribution. Nevertheless an important result that emerges from the lattice calculations is 
that for sufficiently strong coupling such singular configurations are suppressed by quantum fluctu¬ 
ations and thus by the nature of the measure, which imposes nontrivial constraints coming from the 
generalized triangle inequalities. The lattice results suggest therefore that the conformal instability 
is entirely cured for sufficiently strong coupling. It is characteristic of first order transitions that 
the free energy develops an infinitely sharp delta-function singularity at ku, with the metastable 
branch developing no nonanalytic contribution at ku- Indeed it is well known from the theory of 
first order transitions that tunneling effects will lead to a purely imaginary contribution to the free 
energy, with an essential singularity for k > ku US). In the following we shall therefore clearly 
distinguish the instability point ku from the true critical point at kc- Consequently the nonanalytic 
behavior of the free energy (and its derivatives which include, for example, the average curvature) 
has to be obtained by analytic continuation of the Euclidean theory into the metastable branch. 
This procedure is then formally equivalent to the construction of the continuum theory exclusively 
from its strong coupling (small k or large G) expansion, for example starting from 

OO 

ZUk) = Y^ank^ , (64) 

n=0 

OO 

n{k) = > ( 65 ) 

n=0 

OO 

XTlik) = . ( 66 ) 

n=0 

Given a large enough number of terms in this expansion, the nonanalytic behavior in the vicinity of 
the true critical point at kc can then be determined unambiguously, using for example differential 
or Fade approximants mm for suitable combinations which are expected to be meromorphic in 
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the vicinity of the true critical point. In the present case, instead of the analytic strong coupling 
expansion, one makes use of a set of (in principle, arbitrarily) accurate data points to which the 
expected functional form can be fitted. What is assumed here then is the kind of regularity which 
is always assumed in extrapolating finite series to the boundary of their radius of convergence. 
Ultimately it should be kept in mind though that one is really interested in the pseudo-Riemannian 
case, and not the Euclidean one for which such an instability due to the conformal mode is, as 
stated before, to be expected. Indeed had such an instability not occurred one might wonder if the 
resulting theory still had any relationship to the original continuum theory: for the lattice theory 
one expects such an instability to develop at some point, since the continuum theory is known to 
be unstable for weak enough coupling. In conclusion, in the following only data for k < ku will be 
considered; in fact to add a margin of safety only k < 0.051 will be considered throughout the rest 
of the paper. 

To extract the critical exponent 6, one fits the computed values for the average curvature to 
the form of Eq. (H^ . It would seem unreasonable to expect that the computed values for TZ are 
accurately described by this function even for small k, away from the critical point at kc- Instead, 
the data is fitted to the above functional form for either k > 0.02 or k > 0.03. Then the difference 
in the fit parameters can be used as one more measure for the error. In addition, it is possible to 
include a subleading correction of the form 

7^(fe) ~ -An[k,-k + B {kc-kfY > (67) 

fc— 

and use the results to further constraint the uncertainties in the amplitude kc and the exponent 
5 = Au — 1. Using this set of procedures for TZ{k) one obtains on a lattice with sites the following 
set of estimates 


II 

kc = 0.07025(20) 

= 0.357(8) 

(68) 

II 

00 

kc = 0.05811(27) 

1/ = 0.308(16) 

(69) 

L = 16 

kc = 0.06134(11) 

ly = 0.322(6) 

(70) 

CO 

II 

kc = 0.06094(10) 

z/ = 0.320(6) . 

(71) 


Then using the same set of procedures for |7^(A:)|^ (which assumes = 1/3 exactly) one obtains on 
the same lattices 

L = 4 kc = 0.06485(20) (72) 

L = 8 kc = 0.06337(27) (73) 
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Figure 4: Average local curvature TZ{k) as defined in Eq. (fHl) . computed on a lattice with 
16^ = 65,536 sites. Statistical errors (~ 0(10“^)) are much smaller than the size of the symbols. 
The continuous line represents a fit of the form A {kc — k)^ for k > 0.02, with exponent 6 = Av — 1 


L = 16 kc = 0.06377(11) 

(74) 

L = 32 kc = 0.06387(9) . 

(75) 


This last result is presumably the most accurate one, since it is derived form the largest lattice, 
with the highest statistics and the smallest errors on the individual data points. All of these results 
are displayed in Figures 4 to 8, and indicate that the exponent u (and therefore 6) is indeed very 
close to 1/3. Specihcally, Figures 6,7 and 8 show a graph of the average curvature Tl{k) raised to 
the third power; one would expect to get a straight line close to the critical point if the exponent for 
TZ{k) is exactly 1/3. The numerical results indeed support such an assumption, and the linearity 
of the results close to kc is quite striking. The computed data is quite close to a straight line over 
a wide range of k values, providing further support for the assumption of an algebraic singularity 
for Tl(k) itself, with exponent close to 1/3. This last value can be compared to the old estimate 
computed in [iS], i/ ss 0.33. 
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Figure 5: Average local curvature TZ{k) as defined in Eq. (fT4|l . computed on a large lattice with 
32“^ = 1,048,576 sites. Note the change in horizontal scale as compared to the previous figure. 
Statistical errors (~ 0(10“^)) are much smaller than the size of the symbols. The continuous line 
represents a fit of the form A (kc — kY for k > 0.04, with exponent <5 = 4z^ — 1 . 
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Figure 6: Average local curvature TZ{k) on the 16^ lattice, raised to the third power. If5 = z/ = l/3 
exactly, then all the data should fall on a straight line close to kc- The continuous line here represents 
a linear fit of the form A {kc — k) for k > 0.02. Deviations from linearity of the transformed data 
are rather small. 
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Figure 7: Average local curvature TZ{k) on the larger 32^ (1,048,576 sites) lattice, raised to the 
third power. Again \i 5 = v = l/2> exactly, then the points should all fall on a straight line. The 
continuous line represents a linear fit of the form A [kc — k). Deviations from linearity of the 
transformed data are rather small. 
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Figure 8: Volume dependence of the average local curvature \R,{k)\^ on lattices with 4*^, 8^, 16^ 
and 32^ sites. Again, if <5 = i/ = 1/3 exactly, then the data should all fall on a straight line close to 
kc- The continuous line represents a linear fit of the form A {kc — k). The size dependence becomes 
rather small on the larger lattices, unless one moves very close to the critical point. 


30 




8 Curvature Fluctutations 


Figures 9 to 12 show the average curvature fluctuation xnik) defined in Eq. (fTGll . At the critical 
point the curvature fluctuation is expected to diverge, by definition. As in the case of the average 
local curvature TZ{k) analyzed previously, one can extract the critical exponent 5 and kc by fitting 
the computed values for the curvature fluctuation to the form given in Eq. (1451) . And, as for the 
average curvature itself, it would seem unreasonable to expect that the computed values for xn{k) 
are accurately described by this function even for small k, away from the critical point. Instead 
the data has been fitted to the above functional form either for k > 0.02 or for k > 0.03, and the 
difference in the fit parameters is then used as a measure for the error. In addition one can include 
here a subleading correction as well, of the form 

Xn{k) -A^^[k,-k + B{k,-kf]-^^-''> , (76) 

K—ykc 

and use the results to further constraint the errors on the amplitude kc and the exponent 

5 = Ah' — 1. 

One finds that the values for 6 and kc obtained in this fashion are consistent with the ones 
obtained from the average curvature TZ{k), but here with somewhat larger errors, since fluctuations 
are notoriously more difficult to compute accurately than local averages, and require therefore 
significantly higher statistics. Using these procedures one obtains on the largest lattices with 16^ 
and 32^ sites 

kc = 0.05383(102) i/ = 0.350(56) . (77) 

Alternatively, one can use for XTzik) the best estimate for kc obtained earlier from the average 
curvature. This then gives 

1/ = 0.321(12) , (78) 

which is closer to the value obtained from TZ{k). 

Figures 11 and 12 show the inverse curvature fluctuation xizik) on the 16^ and 32'^-site lattices, 
raised to power 3/2. One would expect to get a straight line close to the critical point if the 
exponent for X 7 l{k) is exactly —2/3. The computed data is more or less consistent with a linear 
behavior for k > 0.03, providing further support for an algebraic singularity for XTlik) itself, with 
exponent close to —2/3. Using this last procedure one hnds on the largest (16^ and 32^) lattices 
the improved estimate for the critical point 

kc = 0.06369(84) , (79) 
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Figure 9: Curvature fluctuation XTl{^) on lattices with 16^ = 65, 536 sites. The continuous line 
represents a fit of the form XTlik) = A {kc — for k > 0.02. 
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Figure 10: Curvature fluctuation XTz{k) on the 32^ lattice with 1,048,576 sites. Note the change 
in scale from the previous figure. The line shown is a best fit of the form XTzik) = ^ (kc — k)~^^~^'^ 
for k > 0.04. 
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Figure 11: Inverse curvature fluctuation raised to the power 3/2, on the 16^ lattice; note that the 
data is scaled by a factor of xlOO. The straight line represents a linear fit of the form A (kc — k). 
The location of the critical point in k is consistent with the estimate obtained from the average 
curvature, but with a somewhat larger error. 


which is consistent with the value obtained earlier from TZ^ (see Figures 6 to 8 and related discus¬ 
sion), and suggests again that the exponent v must be rather close to 1/3. 

In order to check the consistency of the results so far, it is possible to analyze the previous 
calculations in a different way. From the definition of the average curvature TZ and curvature 
fluctuation [Eqs. (fT4]l and (fTH)) ]. and the fact that they are both proportional to derivatives of the 
free energy F with respect to k [Eqs. (fT7|l and (fTSl) ]. one notices that their ratio is given by 

2(^^) Xn{k) 


n{k) 




(80) 


The assumption of an algebraic singularity in k for TZ and XTl (Eqs. (I44p and (|45ll l then implies 
that the logarithmic derivative as defined above has a simple pole at kc, with residue 5 = Av — 1 


2(/^) Xn{k) ^ 5 

7Z{k) k^kc k — kc 


(81) 


and the critical amplitudes dropping out entirely for this particular ratio. Eigures 13 and 14 show 
the results for the logarithmic derivative of the average curvature TZ{k), obtained from the data 
shown earlier in Eigures 4 to 12. Using this method on the largest 16“^ and 32^ lattices one finds 


kc = 0.06338(55) v = 0.3356(84) 


(82) 


Note that for the quantity in Eq. (l8T]l only two parameters are fitted, as opposed to three earlier. 
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Figure 12: Inverse curvature fluctuation raised to the power 3/2, on the 32^ lattice; note that the 
data is scaled by a factor of xlOO. The straight line represents a linear fit of the form A (kc — k). 
The location of the critical point in k is consistent with the estimate obtained from the average 
curvature on the same size lattice, here with a larger uncertainty. 

which leads to a slightly improved accuracy. It is encouraging that the above estimates are in good 
agreement with the values obtained previously using the other methods. 

As a further check, it is possible to look at the behavior of quantities when compared directly 
to the average local curvature. Figure 15 shows a plot of the curvature fluctuation xui^) versus 
the curvature Tl{k) (as opposed to k). If the average local curvature approaches zero at the critical 
point (where curvature fluctuation diverges), then one would expect these curvature fluctuations 
to diverge precisely at 7^ = 0. One has from Eqs. and ([15]) 

X7^(7^) ~ A ~ ^ . (83) 

k^kc 

An advantage of this particular combination is that it does not require the knowledge of kc in order 
to estimate ly. Consequently only two parameters are fitted, the overall amplitude and the exponent 
in Eq. ([83]) . Using this method one finds, assuming that the fluctuations diverge at 7^ = 0, 

1/ = 0.3322(71) , (84) 

which is rather consistent with previous estimates. Again the error on ly can be obtained, for 
example, by reverting to more elaborate fits of the type 

y7^(7^) ~ A I 77 + 577^ . (85) 
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Figure 13: Inverse of the logarithmic derivative of the average curvature TZ{k) [defined in Eq. ()8ip ] 
on the 16^ lattice with with 65, 536 sites. The straight line represents a best fit of the form A (kc—k) 
for k > 0.02. The location of the critical point in k is consistent with the earlier estimate coming 
from the average curvature TZ{k) and its fluctuation XTl{k)- From the slope of the line one then 
computes directly the exponent v. 
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Figure 14: Inverse of the logarithmic derivative of the average curvature TZ{k) [defined in Eq. (I8ip ] 
on the 32^ lattice with with 1,048,576 sites. The straight line represents a best fit of the form 
A (kc — k) for k > 0.04. The location of the critical point in k is consistent with the estimate 
coming from the average curvature TZ{k) and its fluctuation XTl{^)- From the slope of the line one 
computes directly the exponent u. 
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Figure 15: Inverse curvature fluctuation, versus the average curvature TZ. Points shown 

here are for the largest lattices. For i/ = 1/3 exactly, y/X'R is expected to be linear in TZ for small 
7 ^. 

Note also that for = 1/3 the exponent simplifies to —2, and one obtains the simple result (see 
also Figure 15) 

xnm A\n\-‘^ . (86) 

One concludes that the evidence so far supports a vanishing average local curvature at the critical 
point, where the curvature fluctuation XTl thus the correlation length ^ [in view of Eqs. (j45h 
and (|42p ] diverge. These results also show some degree of consistency in the values for kc obtained 
independently from TZ{k) and xnik) (Figures 4 to 15). 


9 Finite Size Scaling Analysis 


A further consistency check on the values of the critical exponents is provided by a systematic 
finite size scaling (FSS) analysis, o Indeed the numerical results presented in the previous sections 
have been obtained separately for each lattices of different size. It would be highly desirable if all 
those results could be combined into a single large dataset which then encompasses all the different 


lattice sizes, with consequently a much higher statistical significance. 

A comprehensive review article can be found in the second of | 52l I53 | ; the subject is also covered in numerous 
books on statistical field theory [1511,35] . A systematic field-theoretic derivation of finite-size scaling based on the 
renormalization group is given in |54| . 
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Quite in general, the FSS scaling form for a quantity O diverging like t in the infinite volume 
limit is 


0(L,t) = 


fo 


L 




(87) 


with L the linear size of the system, t the reduced temperature or distance from the critical point, 
fo a smooth scaling function, ^( 00 , t) the infinite volume correlation length and w a correction 
to scaling exponent; but for sufficiently large volumes the correction to scaling term involving uj 
can be safely neglected. In the gravity case one has t ~ \kc — k\, and 0{L,t) is some physical 
average such as the local curvature Tl{k) or its fluctuation XTz{k), with the linear size of the system 
L ~< 1/ . General properties of the scaling function foiu) include the fact that it is expected 

to show a peak if the finite volume value for O is peaked, it is analytic at x = 0 since no singularity 
can develop in a finite volume, and fo{y) ~ for large y for a quantity O which diverges as 

the infinite volume limit. 

The expression in Eq. (|87p is only useful when the infinite-volume correlation length is accu¬ 
rately known. Nevertheless close to the critical point one can use ^ ~ and then deduce from it 
the equivalent scaling from 


0{L,t) = \fo{Le) + 0{L-‘^) 


( 88 ) 


which relies on a knowledge of t, and thus of the critical point, instead. 

The hnite size scaling behavior of the average local curvature, as defined in Eqs. m and (d 
will be discussed next. If scaling involving k and L holds according to Eq. (I88l) . with xq = 1 — 4zz 
the scaling dimension for the curvature, then all points for different k’s and L’s should lie on the 
same universal curve. Erom Eq. ()88p . with t ^ kc — k and xq = —S = 1 — 4z^, one has 


n{k,L) 


^-(4-1/1/) 


n[ikc 


k) + 0{L-^) 


(89) 


where again w > 0 is a correction-to-scaling exponent. The above argument then suggests that the 
quantities 

/ r \ 4-1/1/ 

n{k,L)-L^-^^^ ~ (90) 

should all he on a single universal curve when displayed as a function of the scaling variable 

X ^ (fce-fc)Li/" ~ . (91) 

Eigure 16 shows a graph of the scaled curvature TZ{k) for different values of L = 4,8,16, 32, 

versus the scaled coupling {kc — . The data does indeed support such scaling behavior, and 
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Figure 16; Finite size scaling behavior of the scaled curvature TZ{k, L) ■ versus the scaled 

coupling {kc — k) ■ . Here L = 4,8,16,32 for the lattice with sites. Statistical errors are 

comparable to the size of the dots. The continuous line represents a best fit to a scaling function of 
the form a+bx^^, and finite size scaling predicts that all points should lie on the same universal curve. 
The continuous line corresponds to a critical point kc = 0.06388(32) and exponent = 0.3334(4). 


one finds a best fit for 

k^ = 0.06388(32) z/ = 0.3334(4) 


(92) 


Note that the value for kc found here is in good agreement with the value given earlier in Eq. (I75p . 
Thus so far the finite size scaling analysis lead to values for kc and v which are in good agreement 
with what was obtained before, and provides one more stringent test on the value for v, which 
appears to be again consistent, within errors, with v = 1/3. 

The hnite size scaling properties of the curvature fluctuation, defined in Eqs. ()15jl and ()16ji will 
be discussed next. Again, if scaling involving k and L holds according to Eq. ([ 88 ]) . with t ^ kc — k 
and xo = l — (5 = 2 — 4i/ then all points should lie on the same universal curve. From the general 
form in Eq. (I88p one expects for this particular case 


XTz{k,L) 


^ 2/.-4 


X7e 



k) + 0{L-^) 


(93) 


where a; > 0 again a correction-to-scaling exponent. The above arguments then suggests that the 
quantity 

/ r \ 4-2 /i/ 

Xn{k,L)-L^-^/’^ ~ (^-j (94) 

should give points all lying on a single universal curve when displayed again as a function of 
the scaling variable x in Eq. (fM]l . Figure 17 shows a graph of the scaled curvature fluctuation 
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Figure 17: Finite size scaling behavior of the scaled curvature fluctuation XTl{k, L) ■ versus 

the scaled coupling — k) ■ Here L = 4,8,16, 32 for a lattice with sites. The continuous 

line represents a best fit to a scaling function of the form l/(o + 6x'^), and finite size scaling predicts 
that all points should lie on the same universal curve. The continuous line corresponds to a critical 
point kc = 0.06384(40) and an exponent u = 0.3389(56). 

Xv,{k')ll?‘l'^~^ for different values of L = 4, 8,16,32, versus the scaled variable {kc — k)L^^’^. Using 
this method one finds approximately 

kc = 0.06384(40) i/ = 0.3389(56) . (95) 

Note that the errors in this case are much larger than for the corresponding average curvature 
analysis. Nevertheless the data supports such scaling behavior, and suggests again that v is close 
to 1/3. 

The value of kc itself is expected to have a weak dependence on the linear size of the system 
To ~ . For a hnite system of linear size Lq one anticipates |52l [53] that close to the critical 

point 

kciLo) ~ kc{oo) + c . (96) 

Lq^oo 

This is essentially the expression in Eq. (I42p . with ^ ~ Lq, and then solved for the finite volume 
critical point kc{LQ). Indeed such a weak size dependence is found when comparing kc (as obtained 
from the algebraic singularity fits discussed previously) on different lattice sizes. Figure 18 shows 
the size dependence of the critical coupling kc as obtained on different size lattices. In all three 
cases kc{L(j) is hrst obtained from a fit to the average curvature of the form Tl{k) = A {kc — kY 
as in Eq. (|44p . Due to the few values of L it is not possible at this point to extract an estimate 
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Figure 18: Size dependence of the critical point kc{L) for different lattices, with N = sites and 
L = 4,8,16,32. The line represents a fit kc{L) = kc{oo) + AjL^ + B/L^ and gives the limiting 
estimate kc{oo) = 0.0638615. 

for V from this particular set of data. But since v is close to 1/3, it makes sense to use this value 
in Eq. (|96p . at least as a first approximation. So if one assumes = 1/3 exactly and extracts kc 
from a linear ht to |7^|^, then the variations in k^ for different size lattices are substantially reduced 
(points labeled by smalll circles in Figure 18). This then gives one additional independent estimate 
(which now combines all available lattice sizes, namely L = 4, 8,16,32) 

kc{oo) ~ 0.063862 (97) 

which is in good agreement with the value from the finite size analysis given in Eq. (1921) 

One physicsl qnantity of significant interest is the fundamental gravitational correlation length 
^(/c) itself. It is defined via the exponential decay of physical correlations [such as the ones given in 
Eqs. (f2T)l and ([2^ ] as a function of the geodesic distance between points [see for example Eq. mi 
It also appears as the qnantity of key significance in the scaling argument for the free energy [see 
Eq. (|41ll ]. and is expected to diverge in accordance with Eq. (1421) in the vicinity of the critical 
point at kc- The discussion given in the previous sections pointed to the fact that this qnantity is 
small and of the order of one average lattice spacing (^ ~ /q) in the strong coupling limit (small 
A:), and is expected to increase monotonically towards the critical point at kc in accordance with 
Eq. ()42p . Indeed all the results presented in the previous two sections have been analyzed in terms 
of nniversal scaling properties in accordance with the basic assumption of Eq. gm, and all the 
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results that follow from it. From the results presented so far one concludes that the correlation 
length exponent v defined in Eq. (|42l) is consistent with v = 1/3. 

The next step is to fix the correlation length critical amplitude as well, which is defined in 
Eq. ()42p . The latter is not obtained in an obvious way from any of the results presented so far, 
and requires instead a direct and separate computation of physical correlations at fixed geodesic 
distance, such as the one in Eq. (1271) . These correlations were already computed in [25], and 
additional estimates on the correlation length ^ and can be obtained separately from the size or 
volume dependence of local averages, which is expected to behave, for fixed k ^ kc but close to the 
critical point, as 

^oo(fc) + —(98) 

where here Lq rsj Ei/4 

is a suitably defined linear size of the system. Nevertheless, the overall 
errors for this analysis can be reduced significantly if one assumes = 1/3 exactly (which from the 
previous results on TZ{k) and is known to be a very good approximation), and furthermore if 

one assumes that the correlation length diverges at one and the same critical kc (also determined to 
great accuracy from the previous results for TZ{k) and XTl{k))- The latter set of results was largely 
based on the scaling assumption in Eq. m- Given these simplifying choices one then obtains 

n{k) ~ (fee - k)^’'-^ ~ ~ 1/e, (99) 

fc —>■ fcc 

and also 

xn{k) ~ (kc - ~ ~ e' ■ (100) 

k^kc 

Therefore the two combinations TZ{k) ■ ^{k) and XTzik)/ (k) are expected to approach a constant 
as k —>■ kc- Computing these combinations is so far the most accurate way of determining the 
dependence on k oi ^{k), and in particular for establishing a numerical value for the key amplitude 
in Eq. (Ii2]l . Via this route one finds close to kc that TZ ■ ^ = A^ ■ A-ji ~ 19.57 and xn!^ = 
A^jA^ ~ 2.216, which then gives for the correlation length amplitude in Eq. (|42l) the estimate 
A^ ~ 0.80(3). A plot of the correlation length ^{k) obtained in this way is shown in Figure 19. 
Note that a knowledge of the amplitude A^ then gives immediately, by the renormalization group 
equations in Eqs. (|59p . (|61l) and (|63p . the running of G in the vicinity of the nontrivial fixed point 
at Gc- 
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Figure 19: Estimate for the gravitational correlation length ^(k) versus bare coupling k. For a 
correlation length exponent v = lf3 [see Eq. (H2|) ]. is expected to be linear in k close to 

the critical point kc- 

10 Summary of Results 


Table I summarizes the results obtained for the critical point kc = l/SvrGc and for the universal 
critical exponent u obtained so far using a variety of observables and methods. In view of the 
detailed discussion of the previous section one finds from the best data so far (the one with the 
smallest statistical uncertainties, and the least systematic effects) 

/cc = 0.063862(18) = 0.334(4) , (101) 

which is consistent with the conjecture that v = 1/3 exactly for pure quantum gravity in four 
dimensions. In turn this gives for the bare coupling G at the critical point 

Gc = T-V = 0.623042(25) . (102) 

OTT rVQ 

In previous work [38] the following estimates were given 

kc = 0.0636(11) 1/ = 0.33(1) , (103) 

which have been refined in view of the higher statistics and larger lattices which are part of the 
current study. 
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Observables used to compute kc and u 

Critical Point kc 

Universal Exponent u 

Average Curvature TZ vs. k 

0.06336(28) 

0.331(4) 

Average Curvature vs. k 

0.06367(29) 

0.332(2) 

Average Curvature 7^^ vs. k 

0.06407(24) 

- 

Curvature Fluctuation xn vs. k 

0.05383(102) 

0.350(56) 

Curvature Fluctutation xn vs. k 

- 

0.321(12) 

Curvature Fluctuation Xpi^"^ vs. k 

0.06369(84) 

- 

Logarithmic Derivative 2{1^)x'R.I'TZ vs. k 

0.06338(56) 

0.336(8) 

Curvature Fluctuation xn vs. TZ 

- 

0.332(7) 

TZ{k, L) Finite Size Scaling 

0.06388(11) 

0.333(2) 

XTz{k,L) Finite Size Scaling 

0.06384(18) 

0.339(6) 

Size Dependence of the Critical Point kc{L) 

0.063862(30) 

- 


TABLE I. Summary of results for the critical point kc and the universal gravitational critical exponent v, 

as obtained from the largest lattices studies so far. 


Table II provides a comparison between the best lattice estimate given in Eq. (|IUip and the 
value of the universal exponent v from other approaches. These include the calculation of in the 
2 + e expansion for gravity [22] carried out to two loop order [23| , a calculation of the same using a 
truncated renormalization group approach in four dimensions [551 ES] ; including some recent more 
refined estimates [57|. Further references include a simple geometric argument based on geometric 
features of the graviton vacuum polarization cloud, which gives u = l/(d — 1) for large d |33j . 
and the rough estimate for v from the lowest nontrivial order strong coupling expansion for the 
gravitational Wilson loop [29], which gives v = \l‘l. Finally the result of [46] is mentioned, where it 
was found that a solution to the nonlocal effective field equations of Eq. ()62p for the static isotropic 
metric can only be found provided u = l/{d — \) exactly for d > 4. For a plot of the corresponding 
values for u see Figure 20. 


Table III then gives a similar table, where values for the universal gravitational critical exponent 
V are given in three dimensions (for the Euclidean case) or 2+1 dimensions (for the Lorentzian case). 
Here again it is possible to make a direct comparison between several approaches, namely the lattice 
[58[ Wf\ . the two-loop 2 + e expansion of [23] but with now e = 1, the Einstein-Hilbert truncated 
renormalization group approach [55] [5(3 E] and the large d estimate of m- The corresponding 
values for u are shown in Figure 20. 
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Method used to compute v in d=4 

Universal Exponent u 

Euclidean Lattice Quantum Gravity (this work) 

= 2.997(9) 

Perturbative 2 + e expansion to one loop |22] 

= 2 

Perturbative 2 + e expansion to two loops [23] 

= 22/5 = 4.40 

Einstein-Hilbert RG truncation [56] 

PS 2.80 

Recent improved Einstein-Hilbert RG truncation [57| 

^ 3.0 

Geometric argument [33] Pvacpoi{r) ~ r'^~^ 

= d — 1 = 3 

Lowest order strong coupling (large G) expansion [29| 

= 2 

Nonlocal field equations with G(n) for the static metric |l6| 

= d — 1 for d > 4 


TABLE II. A comparison of estimates for the fundamental scaling exponent v, based on a variety of 
different analytical and numerical methods. These include the 2 + e expansion for pure gravity carried out 
at one and two loops \2^ . an estimate for the leading exponent in a truncated renormalization group 
expansion J561 157| /, a simple geometric argument based on the geometric features of the quantum vacuum 
polarization cloud for gravity, and finally the only value allowed by a consistent solution to the nonlocal 
field equation with a G'(n) for the static isotropic metric. 


Method used to compute in d = 3 

Universal Exponent u 

Euclidean Lattice Quantum Gravity [58| 

= 1.72(5) 

Exact solution of Lorentzian Gravity (Wheeler-DeWitt Eq.) in 2-|-l dim. |47j 

= 11/6 = 1.8333 

Perturbative 2 -|- e expansion to one loop [22| 

= 1 

Perturbative 2 -|- e expansion to two loops [23] 

= 8/5 = 1.6 

Einstein-Hilbert RG truncation [56] 

K. 1.33 

Large d geometric argument [33] Pvacpoiir) ~ 

= d-1 = 2 


TABLE III. A comparison of various estimates for the fundamental scaling exponent v in 2 + 1 
dimensions, based on a variety of different analytical and numerical methods. Included are the 2 + e 
expansion for pure gravity carried out at one and two loops 1221 \2Stj . an estimate for the leading exponent 
in a truncated renormalization group expansion and a simple geometric argument based on general 
features of the quantum vacuum polarization cloud for gravity. 
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Figure 20: Universal scaling exponent v determining the running of G [see Eqs. (1591) and (1611) ] as 
a function of spacetime dimension d. Shown are the results in 2 + 1 dimensions obtained from the 
exact solution of the lattice Wheeler-DeWitt equation m, the numerical result in four dimensions 
(this work and [H] ), the 2 + e expansion result to one [22] and two loops [23] , and the large d result 
~ d — 1 [33] . For actual numerical values see Tables II and III. 

11 Implications from a Gravitational Exponent z/ = 1/3 


In this section the consequences of having a definite value for the critical point kc, as well as a 
value for the universal critical exponent v [see Eq. (jlOlj . will be discussed. One notes on the one 
hand that the value for kc and thus Gc now fixes the value for the ultraviolet cutoff a. At the same 
time, the specific value for i/ ~ 1/3 gives predictions for the scaling behavior of local averages, 
gravitational correlations and the running of G. 

Eirst note that the value for the critical point kc given in Eqs. (llOip and (|102l) fixes the lattice 
spacing a, and thus the value for the ultraviolet cutoff 

G ^ Gc = 0.623041 . (104) 

Erom the known laboratory value of Newton’s constant G, Ip = ■sJhGj^= 1.616199(97) x 10“^^cm, 
so that one obtains for the fundamental lattice spacing a = 1.2669 y/G^ = Ip, or 

a = 2.0476 x 10"^^ cm . (105) 

and from it a definite value for the ultraviolet cutoff A ~ 1/a. Eor the average lattice spacing in 
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units of a one finds 

<P>= 1“^ =[2.398{9)af (106) 

so that a and Iq are quite comparable in magnitude [this fact can be traced back to the original 
overall scale choice Aq = 1 in Eqs. ([2]) and ([6]), motivated by Eq. (fTO]) ]. 

For the average local curvature TZ{k) one has from Eq. (l44h . using Eq. (j42h and ly = 1/3, 


< f d‘^x y/g R(x) > 




1/iy—d 




(107) 


< I d^Xy/g> 

The dimensionless amplitude is expected to be 0(1) in lattice units, and is given below in 
Eq. (11111) . This result is based on the fact that the lattice calculations allow one to also extract 
various amplitude coefficients. For the dimensionless curvature amplitude defined in Eq. (I44|) one 
finds 

An = 24.46(9) , (108) 

and for the dimensionless curvature fluctuation amplitide defined in Eq. (I45p 

Av-l 


Ay^ — 


■An = 1.418(6) 


</2> 

Combined with the dimensionless correlation length amplitude defined in Eq. 


(109) 


= 0.80(3) 


( 110 ) 


one finds for the amplitude in Eq. (|107p 

= 3-40(13) . 

For the curvature fluctuation Xnik) one has from Eqs. (fTCIl . (1451) and p00|) 

< (J cAx y/g R)"^ > — < f d^x y/gR 
< f d^Xy/g> 

with dimensionless amplitude 




( 111 ) 


( 112 ) 


jsi = hi. 
^ - A} 


Av — 1 An 


<P> 




= 2.22(9) 


(113) 


These results in turn provide useful information for the curvature correlation function at fixed 


geodesic distance of Eqs. (l24j) and ([2^ . From Eq. 
2n = 2(fi-1/z^) = 2(4-3) = 2 


one has for the power appearing in Eq. 


In weak field perturbation theory one finds [26] < y/gR{x)y/gR{y) >c ~ < d^h{x)d^h{y) > ~ l/\x — y|®, so the 
result here is quite different. If one defines an anomalous dimension g for the graviton propagator in momentum space, 
< /i /i >~ one finds ?7 = d—2 — 2/i/orr7 = —4in four dimensions for i/ = 1/3, which deviates significantly from 

the gaussian or perturbative value. Such a large deviation is already observed in the 2 + e expansion [see Eq. mi 
and is not peculiar to lattice quantum gravity. In gravity such a possibility was already discussed some time ago in 

m- 
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One then obtains for the curvature-curvature correlation function at “short distances” r <C ^ 
and for z/ = 1/3 the remarkably simple result 

1 ^0 


< ^/g Rix) y/g Riv) 5(|a: -y\-d) >c 


(114) 


d2d-2/y a2 (p 

Note that in the last term the correct dimensions have been restored, by inserting suitable powers 
of the lattice spacing a. It is instructive to compare the above result to the expression for the local 
average curvature, Eq. (|107p : note in particular that both expressions still contain explicitly the 
size of the microscopic parallel transport loop ~ a ~ Zp. Here the dimensionless amplitude Aq is 
related to the amplitude in Eq. (jll3p because of Eq. (f47p . and one finds 


. _ _ 1 1 
° “ 27r2 “ 27r2 ' 3 


1 


_ 

<P> ' Aj 


= [0.335(20)]^ 


(115) 


so that the dimensionless correlation function normalization constant is = y/A^ = 0.335(20). 
As expected, all of these amplitudes are close to 0(1) in units of the ultraviolet cutoff (fundamental 
lattice spacing) a. 

The exponent p = 1/3 and the amplitude A^ now determine the running of G with scale, see 
Eq. ([5^ . and one obtains 


G{q^) = Gc 


1 + Co 


m 


2\ 3/2 


+ 0( 


m 


(116) 


with reference scale m = 1/^. The coefficient cq [see Eg. (1601) ] determines the amplitude of the 
quantum correction, and it is given by 

CO = SvrOeAj/" = Al/kc , (117) 

with = 0.80(3) from the numerical solution [see Eq. (IllOp : for the definition of the amplitude 
A^ see Eq. (|42l) ]. and also = 1/3 and kc from Eq. (|101l) . This then gives for the dimensionless 
amplitude of the leading quantum correction in Eq. (I116h co ~ 8.02. This last result can then be 
translated directly into a covariant G'(n) [see Eq. (l6Tp with u = 1/3], 

1 


G'(n) = G, 


1 + Co 


\ 3/2 

d +... 


( 118 ) 


_^2q2^ 

The latter forms the basis for a set of nonlocal effective field equations [see Eq. (I62P ]. and gives the 
running of G{r) for the specific choice of the static isotropic metric [see Eq. (|63p ]. 

In principle it is also possible to estimate the next (sub-leading) correction to the leading running 
of G in Eq. (11161) . given the knowledge of the subleading corrections to ^{k) or m{k) in Eqs. 
or ()116p . If one has 

^-\G) = a-iA„(G - G,f [1 - hu\G - G,) + 0((G - G,f)] , 
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with a sub-leading correction of amplitude b, then for the running of G one obtains to this order 


G{q^) 

Gc 


1 + Co 


V 


+ Cl 



l/v 


-I- . . . , 


( 120 ) 


with Cl = hvj{J^rn Gc) ~ 2.87, given that h ps 0.215 and Am = {kdGcYjA^, and A^ given in 
Eq. (IllOp . The domain of validity for the above expression isgS>m = l/^orr<C^; the strong 
infrared divergence at g ~ 0 is largely an artifact of the current expansion, and can be regulated 
either by cutting off the momentum integrations at g ~ m = l/(^, or by the replacement on the 
r.h.s. —)• + m?. 

Furthermore, the previous results show clearly that the reference scale for the running of G 
is set by the correlation length which by Eqs. (H6|) . (IM|) and Eqs. (I107p appears to be directly 
related to curvature. In particular the form of the running of G with scale suggests that no 
detectable corrections to classical gravity should arise until either a) the scale r approaches the 
very large (cosmological) scale or b) until one reaches extremely short distances comparable to 
the Planck length r ~ /p, at which point higher derivative terms, light matter corrections and string 
contributions come into play. In other words, the results of Eqs. ([63]) (I116p or ()118|) imply that 
classical gravity is largely recovered on atomic, laboratory, solar and even galactic scales, as long 
as the relevant distances satisfy r <C 

Therefore one crucial ingredient needed in pinning down the magnitude of the quantum cor¬ 
rection for G{(Y‘) in Eqs. (11161) or (jllSp is the actual value of the nonperturbative reference scale 
It was argued in [29| that, in analogy to ordinary gauge theories, the gravitational Wilson loop 
provides precisely such an insight. The main points of the argument are rather simple and can be 
reproduced in a few lines. In complete analogy to the gauge theory case, these arguments basi¬ 
cally rely on the concept of universality, the existence of a universal correlation length at strong 
coupling, and the use of the Haar invariant measure to integrate over large fluctuations of the 
fundamental local parallel transport matrices. Following [261127| . in [29] the vacuum expectation 
value corresponding to the gravitational Wilson loop was defined as 


< W{G) > = < tr [u{G) Ui U 2 ... Un] > 


( 121 ) 


Here the I/’s are elementary rotation matrices, whose form is determined by the affine connection 
and which therefore describe the parallel transport of vectors around a loop C, see also Eq. (|32l) . 
Here oj^y{G) is a constant unit bivector, characteristic of the overall geometric orientation of the 
loop, giving the normal to the loop. In the continuum the combined rotation matrix U(C') is given 
by the path-ordered {V) exponential of the integral of the affine connection T^^, as in Eq. (1311) . so 
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that the previous expression represents a suitable regularized and discretized lattice form. In 
it was then shown that quite generally in lattice gravity for sufficiently strong coupling one obtains 
universally an area law for near planar loops, 


< W(C) > ~ exp (— Ac!^) 


( 122 ) 


where Ac is the geometric area of the loop. This last result relies on a modified first order formalism 
for the Regge lattice theory |34j . in which the lattice metric degrees of freedom are separated out 
into local Lorentz rotations and tetrads. Moreover, the result of Eq. (jl21l) appears be universal 
since and was shown to hold in all known lattice formulations of quantum gravity in the strong 
coupling regime. In [29] an expression for the correlation length ^ appearing in Eq. (I122p was given 
in the strong coupling limit, where one finds ^ = A/■sJkc\log{k/kc)\ + 0{k^). For k close to kc this 
gives immediately ^ ~ A:\kc — k\~^/'^ and thus, to this order, v = ^ and also A^ = A in Eq. (142]) . 
Nevertheless, the discussion of the previous sections and the numerical solution of the full lattice 
theory suggests that the correct expression for ^ to be used in Eq. (I122p should be the one in 
Eq. (|42p . with = 1/3 [Eq. (|lUip ]. kc given in Eq. (|lUip and amplitude A^ = 0.80(3). 

The next step is to make contact between the above results and a semiclassical description, 
which requires that one connects the nonperturbative result of Eq. (|122l) to a suitable semiclassical 
physical observable. Indeed, by the use of Stokes’s theorem, semiclassically the parallel transport of 
a vector round a very large loop depends on the exponential of a suitably coarse-grained Riemann 
tensor over the loop. In this semiclassical picture one has for the combined rotation matrix U 


UMC) 


exp 


'S{C) 




where A^^ is an area bivector associated with the loop in question. 


Ac^ = k j ■ 


Then the semiclassical procedure gives for the loop in question 


(123) 


(124) 


expJ i / R '-Xa^C 

[ Js{C) 

A similar result is of course well established in non-Abelian gauge theories, and by now regarded as standard 
textbook material [see for example, Peskin and Schroeder, An Introduction to Quantum Field Theory, p. 783, Eq. 
(22.3) [60]]. There ^ represents the gauge field correlation length, defined, for example, from the exponential decay 
of connected Euclidean correlations of two infinitesimal chromo-magnetic loops separated by a given distance ja;]. 
Eollowing |29| . we choose to write here the gravitational result in the same scaling form, involving the invariant 
gravitational correlation length an overall, in principle calculable, multiplicative constant 0(1) in the exponent 
has been set equal to one here. 



W{C) ~ tr UJ{C) 
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Here again uJny{C) is a constant unit bivector, characteristic of the overall geometric orientation 
of the parallel transport loop. By carefully comparing coefficients for the two area terms [29] one 
then concludes that the average large-scale curvature is of order at least in the strong 

coupling limit considered in the cited references. Since the scaled cosmological constant can be 
viewed as a measure of the intrinsic curvature of the vacuum, the above argument then gives a 
positive cosmological constant for this phase, corresponding to a manifold which behaves as de Sitter 
(A > 0) on large scales [29]. These arguments then lead to the suggestion that the macroscopic 
(semiclassical) average curvature is related to ^ by 

{ large scales ~ +i/e , ( 126 ) 


at least in the strong coupling (large G) limit. It is important to note here that the result of 
Eq. (|126p applies to parallel transport loops whose linear size rc is much larger than the cutoff, 
T~c ^ Ip, CL] neverthless in this limit the answer for the macroscopic curvature in Eq. (|126ll becomes 
independent of the loop size or its area [29|. Furthermore, these arguments lead, via the classical 

field equations, to the identification of with the observed (scaled) cosmological constant \ohs, 
18 


1 


■j Aofts ~ -j- 


(127) 


3 -oos - . ^2 

In this picture the latter is then regarded as the quantum gravitational condensate, a measure of 
the vacuum energy, and thus of the intrinsic curvature of the vacuum. Then a suitable effective 
action, describing the residual effects of quantum gravity on very large distance scales, is of the 
form 


Ies\9p.i'\ — 1 Gvr G(/i) j ^ ^ y/d ^ 2 ) Imatter [dpu, ■ ■ ■] , (^28) 

with G(/i) a very slowly varying (on macroscopic scales) Newton’s constant, in accordance with 
Eqs. (I116P or (|118p . 

Note that the above results in many ways parallels what is found in non-Abelian gauge theories, 

where for example one has for the color condensate < >~ Furthermore, this last result 

can also be obtained largely from purely dimensional grounds, once the existence of a fundamental 

correlation length which for QCD is given by the inverse of the mass of the lowest spin zero 

glueball, is established. Accordingly, for gravity too one would expect, again simply on the basis of 

^®up to a constant of proportionality, expected to be of order unity. 

Note that, quite generally and independent of the lattice results, it seems rather difficult to implement a weakly 
running cosmological constant, if general covariance is to be maintained at the level of the effective field equations. 
If the running of A is formulated via a A(n) then because of — 0 one also has = 0, which makes it 

nearly impossible to have a nontrivial A(n) |24| . 
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Figure 21: Running gravitational coupling G{r) versus r, obtained from the G{k) in Eq. (1591) by 
setting q ^ 1/r, with exponent u = l/2> and amplitude oq — 8.02(55). The lattice quantum gravity 
calculations done so far suggest roughly a 5% effect on scales of 0.187 x 4890Mpc 910Mpc, and 

a 10 % effect on scales of 0.238 x 4890Mpc llOOMpc. 

dimensional arguments, that the large scale curvature (the graviton condensate) should be related 
to the fundamental correlation length by (i?) ~ 1/^^, as in Eq. (|46l) . 

The considerations presented so far can to some extent finally provide a quantitative handle on 
the physical magnitude of the nonperturbative scale Erom the observed value of the cosmological 
constant one obtains an estimate for the absolute magnitude of the scale ^ 

^ ~ pa 4890Mpc . (129) 

Irrespective of the specific value of this would indicate that generally the recovery of classical GR 
results only happens for distance scales much smaller than the correlation length In particular, 
the Newtonian potential acquires a tiny quantum correction from the running of G{r), 

V{r) = -G(r).^^^ , (130) 

r 

with G{r) given, for the static isotropic solution, in Eq. (|63p . and for which quantum effects become 
quite negligible on distance scales r Eigure 21 shows the expected qualitative behavior for 

the running G{k) over scales slightly smaller or comparable to with the main uncertainty arising 
from estimating the physical magnitude of ^ itself [Eq. (I129p ]. Specifically, from Eq. (I59p the lattice 
prediction at this point is for roughly a 5% effect on scales of 0.184 x 4890Mpc ~ 900Mpc, and a 
10 % effect on scales of 0.232 x 4890Mpc 1130Mpc. 
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The above results also suggest that the curvature on very small scales behaves rather differently 
from the curvature on very large scales, due to the quantum fluctuations eventually averaging out. 
Indeed when comparing the result of Eqs. (|46p and (j99p to the one in Eq. ()126p one is lead to 
conclude that the following change has to take place when going from small (linear size ~ Ip) to 
large (linear size S> Ip) parallel transport loops 


{R) 


small scales 


ipC 


{ ^)large scales 


(131) 


An intuitive way of understanding the above result is that on small scales the strong local fluctua¬ 
tions in the metric/geometry lead to large values for the average rotation of a parallel-transported 
vector. But then on larger scales these short distance fluctuations tend to average out, and the 
combined overall rotation is much smaller, by a factor 0{lp/^), 


Zr = 


II. 

1 


(132) 


The above quantity should then be regarded as an essential and necessary “renormalization con¬ 
stant” when comparing curvature on different length scales, and specifically when going from very 
small (size ~ Ip) to large (size ^ Ip) parallel transport loops. See also the earlier discussion 
preceding Eq. ([351) . about the issue of comparing correlations of large loops versus correlations of 
small (infinitesimal) loops. 

To conclude this section, one can raise the legitimate concern of how these results are changed 
by quantum fluctuations of various matter fields; so far all the results presented here apply to 
pure gravity without any matter fields. Therefore here and in the rest of the paper what has 
been followed is the quenched approximation, wherein gravitational loop effects (perturbative and 
nonperturbative) are fully accounted for, but matter loop corrections are entirely neglected. In 
the presence of matter fields coupled to gravity (scalars, fermions, vector bosons, spin-3/2 fields 
etc.) one would expect, for example, the value for v to change due to vacuum polarization loops 
containing these fields. A number of arguments can be given though for why these effects should 
not be too dramatic, unless the number of light matter fields is rather large. Firstly one notices 
that to leading order in the 2 -|- e expansion the exponent u only depends on the dimensionality of 
spacetime, irrespective of the number of matter fields and of their type [22], i/ ~ l/(d — 2). Also, 
one can show that in the 2 -|- e expansion for gravity [221123] matter loop corrections appear later in 
the form of factors oc (25 — c) in the /3-function, where c is the central charge corresponding to those 
massless matter fields. In four dimensions the correction is thought to be even smaller oc (48 — c) 
m- Thus unless c is rather large, the matter contribution is quite small even at next-to-leading 
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order in the 2 + e expansion |22l I23j . In addition, in the case of lattice gravity the effects of a 
single light scalar held are so small that they are barely detectable in the numerical evaluations of 
the path integral. In general one would expect signihcant infrared modihcations to gravity coming 
from particles that are either massless or very nearly massless. The evidence so far would therefore 
suggests that the approximation in which vacuum polarization effects of light matter helds are 
entirely neglected should still be useful, at least as a hrst step. 


12 Gravitational Scaling Dimensions and Phenomenology 


The previous section dealt with the fact that one of the main implications of quantum gravity is the 
running of the gravitational constant with scale, in accordance with Eq. ()61l) . There are additional 
consequences which arise from the fact that in general gravitational correlations don’t follow free 
field (gaussian) predictions. One example is the curvature correlation function of Eqs. (f2H) . (I26l) 
and ()48p . for which the final form at this stage is given in Eq. (I114p . This section deals therefore 
with a discussion of the implications of the result given in this last expression for the curvature 
correlation function of Eq. (11141) . specifically with lattice spacing a from Eq. (llOSp and amplitude 
given in Eq. (jllSp . Nji ss 0.335. 

There is clearly a rather substantial difference in scale between the curvature appearing in 
Eq. (j24h and therefore in Eq. (jll4h . and the curvature in Eq. (I134p . In the hrst case the curvature 
involves the parallel transport of vectors around inhnitesimal loops, whose size is determined by 
the ultraviolet cutoff [the lattice spacing, comparable to the Planck length because of G Ri Gc due 
to the slow running of G, with Gc given in Eq. (|102ll ]. In the second case the curvature in question 
refers instead to the semiclassical domain, as described by a set of effective long distance held 
equations, for which the curvature is obtained operationally from the parallel transport of vectors 
around macroscopic loops, of linear size much larger than the Planck length. Such a concern should 
be kept in mind when transitioning from the microscopic result in Eq. (jll4l) to the semiclassical 
result written down below in Eq. (|143l) . 

Eirst consider what can be stated purely at the classical level. One can use the held equations to 
directly relate the local curvature to the local matter mass density. Prom Einstein’s held equations 

R^J.u-\g^J.uR = SttGT^u (133) 
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for a perfect fluid one then obtains for the Ricci scalar, in the limit of negligible pressure, 

R{x) ~ SttG p{x) . (134) 

This last result then allows one to relate local fluctuations in the curvature 5 R{x) to local fluc¬ 
tuations in the matter density 6p{x), which could potentially provide a useful connection to the 
quantum result for the correlation function in Eq. (|26p . Of course, in the Newtonian limit the 
above result simplifies to Poisson’s equation 

A/ioo(x, t) = 8TrGp{x,t) , (135) 


where /iqo = 2(/> and p are the macroscopic gravitational field and the macroscopic mass density, 
respectively. 

In the cosmology literature it is customary to describe matter density fluctuations in terms of 
the density contrast correlation function 


G{r) = < 5p{r) 5p{Q) > . (136) 

The latter is related to its Fourier transform P{q) by 

G{r) = ^ tdqq^P{qf-^, (137) 

2vr^ qr 

and the above expression has to contain both an infrared regulator (p) and an ultraviolet cutoff 
(A), to make sure the integral converges. If the power spectrum P{q) is described by a simple 
power law of the form 

m = ^ , ( 138 ) 

(where n = — s is commonly referred to as the spectral index), then one finds in the scaling regime 
l//ri^rS> 1/A for the density contrast correlation function in real space 

G{r) = CsaopA^~'' > (139) 

with Cs = r(2 — s) sin(7rs/2) / 27r^. Not unexpectedly, the answer appears to be quite sensitive to 
the choice for the ultraviolet and infrared cutoffs. For the specific value s = 1 one has P{q) = ao/q, 
and this then gives 

Gir) = {s = l) , (140) 

the cosmology literature the (dimensionless) galaxy matter density two-point function is usually referred to as 
^(r), but here we want to avoid a possible confusion with the gravitational correlation length T 
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which would seem to reproduce the result in Eq. (|114l) . In practice the observational data for such 
matter density correlations is often presented in the simple form 




(141) 


with exponent 7 and scale tq fitted to astrophysical observations. For an exponent 7 close to two, 
one has by comparing Eq. (|140l) to Eq. (I14ip oq = 27 r^/irQ/A, which still requires a choice of 
cutoffs ^ and A; for the most obvious choice here, namely ^ ~ 1/^ and A ~ l//p, one obtains 


«o — 


27r^ Ip 


(142) 


It is rather tempting at this stage to try to connect the observational result of Eq. (jl41h to the 
quantum correlation function in Eq. (|114p . One then expects for the matter density fluctuation 
correlation also a power law decay of the form 

1 1 1 


< 5p(x,t) 5p{y,t') > 


(143) 


|x-y|<€ a 2 (t) |x - y| 

where a{t) here represents the scale factor. Also, this last correlation function can be made dimen¬ 
sionless by suitably dividing it by the square of some average matter density Po ~ Pc = 3Hq/8ttG. 
By comparing coefficients in Eqs. (|114p and (11411) one finds 7 = 2, and for the length scale in 
Eq. (fTiB 


1 


ro = 




(144) 


Stt G po a 

with vF) — 0.335 the dimensionless amplitude for the curvature correlation function of Eq. (I114p . 
and a the lattice spacing given in Eq. (llOSp . 

The preceding argument nevertheless still contains a fundamental flaw, related to the use, at 
this stage in unmodified form, of the curvature correlation function result of Eq. (|114l) . As already 
discussed previously, that form applies to the correlation of infinitesimal (Planck length or cutoff 
size) loops, which would not seem to be appropriate for the macroscopic (or semiclassical) parallel 
transport loops, such as the ones that enter the field equations Eqs. (|133p and (|134p . and which 
thus relate locally the macroscopic 6R{x) to the 6p{x). It would then seem desirable to be able 
to correct for the fact that the parallel transport loops sampled in Eq. (|134l) are much larger than 
the infinitesimal ones sampled in the correlation function in Eq. (I114|) . As in Eqs. (j38|) . (I13ip and 
(jl32h . the transition to macroscopic loops (linear size ^ a) can be affected in Eq. (I114j) by the 
replacement of ^ This then gives for large (macroscopic size ^ a) parallel transport loops 

Ai 


< Vg R{x) ^/g R{y) (5(|x -y\-d) > 


In weak field perturbation theory one finds < p{x) p{y) >c 
result here is quite different. 


d < 5 
< d'^h(x) d^h(y) > 


(145) 

l/\x — j/|®, so again the 
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with the expectation of a comparable amplitude Ai ~ Aq. This last result then leads to the 
following improved estimate for the macroscopic matter density correlation of Eq. (|136ll 


so that comparing to Eq. (11411) one finds for the exponent 7 = 2 , 
improved value 


(146) 


and for the length scale rg the 


which seems more in line with observational data. For the Fourier amplitude og in Eq. ()138ll one 
has now 


ag = 27T^ 


Ip 


2.85 Ipi . 


(148) 


Observed galaxy density correlations give indeed for the exponent in Eq. (11411) a value close to two, 
namely 7 ~ 1.8 ± 0.3 for distances in the O.lMpc to bOMpc range [621163] . and for the length scale 
rg ~ lOMpc. More recent estimates for the exponent 7, going up to distance scales of lOOMpc, 
range between 1.79 and 1.84 [Ml[65lEU|671 |68l[69] . Nevertheless at this point the (perhaps rather 
naive) identification given in Eqs. (I147p and (I148j) . while intriguing, is possibly entirely accidental 
since it bypasses any concerns about the actual physical origin of the galaxy correlation function in 
Eq. (11441) . including the form and evolution of primordial density perturbation, the detailed nature 
of linear relativistic density perturbation theory for a given comoving background etc. 


13 Conclusions 


In this work a number of improved estimates have been presented for gravitational scaling dimen¬ 
sions and amplitudes, obtained from the lattice theory of gravity. Numerical methods combined 
with modern renormalization group arguments and finite size scaling have been shown to provide 
detailed information about rather subtle nonperturbative aspects of the theory. It has been known 
for some time that the Euclidean lattice gravity theory has two phases, only one of which, the 
gravitational anti-screening phase for G > Gc, is physically acceptable. Here we have described in 
some detail the properties of the latter smooth phase, and provided quantitative estimates for the 
critical point, scaling dimensions and the behavior of physical correlations for distances large com¬ 
pared to the lattice cutoff. In many ways the present calculation is still incomplete, in particular 
the gravitational Wilson loop and the correlation between loops has not been studied numerically. 
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and only some general properties have been inferred. Also, more heavy work is needed to accurately 
determine the curvature correlation functions versus distance, and from it the fundamental non- 
perturbative correlation length and various amplitudes connected to it. Furthermore, the derivation 
of a number of results has relied heavily on basic renormalization group scaling, with only a handful 
of explicit checks. Nevertheless, it would seem from the results presented so far that the feasibility 
of these types of calculations should increase signihcantly in the near term due to expected rapid 
advances in hardware and software tools. 

It is encouraging that four different approaches to quantum gravity give rather comparable 
results for the scaling dimensions (see the comparison Tables II and III, as well as Figure 20), 
and therefore suggest a unique underlying renormalization group universality class, associated with 
the quantum version of General Relativity. It is characteristic of the model described here that 
the growth of G with scale is described by a nonperturbative correlation length related to 
the gravitational vacuum condensate, for which a specific quantitative estimate was given earlier. 
More generally, the vacuum condensate picture of quantum gravity presented in this paper makes in 
principle a number of specific and testable predictions, which could either be verified or disproven 
in the near future, as new and increasingly accurate satellite observations become available. The 
main aspects of this picture can be summarized as follows: 

o The vacuum condensate picture of quantum gravity contains from the start a very limited 
number of parameters, and is therefore rather strongly constrained. While it does involve a 
new nonperturbative scale (the gravitational vacuum condensate), it is found that this scale 
simultaneously determines the running of G with scale, the value of the scaled cosmological 
constant, and the long distance behavior of physical invariant correlations. 

o The current theory predicts a slow increase in strength of the gravitational coupling when 
very large, cosmological scales are approached [see Eqs. (|116p and dllSp ]. In this context, the 
observed scaled cosmological constant A acts as a dynamically induced infrared cutoff, similar 
to what happens in non-Abelian gauge theories. In principle, both the universal power and 
amplitude for this infrared growth are calculable with some accuracy from the underlying 
lattice cutoff theory. 

o For a sufficiently large scale ^ (and therefore small A) no observable deviations from classical 
General Relativity are expected on laboratory, solar systems and even galactic scales [see 
Eqs. ([63]) and ()130l) ]. 
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o The calculations presented here give a number of predictions for the behavior of invariant 
curvature correlations as a function of geodesic distance, and specifically the powers and 
amplitudes involved [see Eqs. (jll4|] and (jl45p ]. 

o The lattice theory appears to exclude at this point the possibility of a physically acceptable 
phase with gravitational screening; such a (weak coupling) phase in the lattice theory appears 
to be inherently unstable, presumably as a consequence of the conformal mode, and cannot 
lead to a semiclassical regime for gravity. It leads instead to a pathological degenerate ground 
state describing some sort of branched polymer. Nevertheless for large enough quantum 
fluctuations (large G) the instability is overcome and a new stable phase emerges. 

o In the strong coupling limit (for the Euclidean case) of the lattice theory the effective, long 
distance cosmological constant is positive [29]. In this same regime it seems nearly impossible 
from the lattice theory to get a negative value for this quantity, irrespective of the choice 
of boundary conditions (which in the lattice context play no role in the argument). Also, a 
positive cosmological constant is interpreted here as a genuinely nonperturbative gravitational 
vacuum condensate [see Eqs. (jl26p and (jl27l) ]. 

If the picture presented in this paper is indeed close to correct, then it points to what appears 
to be a deep analogy between the nonperturbative vacuum state of quantum gravity and known 
properties of strongly coupled non-Abelian gauge theories (or what could be called the QCD anal¬ 
ogy). Over time this analogy has been helpful in illustrating properties of quantum gravity, many 
of which are ultimately based on rather basic principles of the renormalization group, connected 
with the scaling properties expected in the vicinity of a nontrivial fixed point. Indeed in QCD there 
exists also a nonperturbative mass parameter m = I/^ (sometimes referred to as the mass gap) 
which is known to be a renormalization group invariant; that such a mass scale can be generated 
dynamically is known to be a highly nontrivial outcome of the renormalization group equations 
for QCD. Eurthermore, there seems to be a fundamental relationship between the nonperturbative 
scale ^ (or inverse renormalized mass) and a nonvanishing vacuum condensate for these theories, 
both for gravity and QCD, 

^ = ( 149 ) 

An additional relevant example that comes to mind is the fermion condensate in gauge theories, 

(V'V’) - ^ , (150) 
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a consequence of confinement and chiral symmetry breaking. Current lattice and phenomenological 
estimates for QCD cluster around ( ^ ) ~ (440MeC)^ and () ~ f290Mel4l^ [70lI7T1172] . 

Modifications to the static potential in gauge theories are best expressed in terms of the running 
coupling constant 0 : 5 (^), whose scale dependence is determined by the celebrated beta function 
of QCD with coupling as = On the one hand, a solution of the renormalization group 

equations give for the running of as{g) 


asifi) = 


dvr 


/3oln/i2/A|^ 


+ ... 


(151) 


On the other hand, the nonperturbative scale ^j^s appears as an integration constant of the 
renormalization group equations, and is therefore - by construction - scale independent. 


^MS = A exp ( - 


/ 


as ( A) 



(152) 


where here A represents the QCD ultraviolet cutoff. The physical value of cannnot be 

fixed from perturbation theory alone and has to be determined instead from experiment, A^^ ~ 
210Mel4. In quantum gravity the corresponding statements are given in Eqs. (|51ll and (1591) . 

Wilson loop correlations play an important role in QCD as they do in quantum gravity. In non- 
Abelian gauge theories a confining potential is found at strong coupling by examining the behavior 
of the Wilson loop, defined for a large closed loop C as 


{W{C)) = (trFexpji^y A^{x)dx^^) , (153) 

with and the fa’s the group generators of SU{N) in the fundamental representation. In 

the pure gauge theory at strong coupling, the leading contribution to the Wilson loop is known to 
follow an area law for sufficiently large loops. The analogous quantity for gravity is the gravitational 
Wilson loop described, for example, in Eqs. (f3Tl) and ()12ip . But in contrast to QCD, in gravity the 
Wilson loop bears no relationship to the static potential [30] (the path ordered line integral of the 
affine connection does not in any way describe a gravitational interaction energy). 

A central role in this view of quantum gravity is played by the gravitational correlation length 
of Eqs. (1271) . ()i 2 ]) and (f52]l . Gauge theories also contain a nonperturbative, dynamically gener¬ 
ated quantity the gauge held correlation length, and it is essentially the same [up to a factor 
0(1)] as the inverse of The same universal quantity ^ also appears in a number of other 

physical observables, including the exponential decay of the Euclidean correlation function for two 
inhnitesimal loop operators separated by a distance |x|, 

Gioop-ioop{x) = {trVex.p^ig I A^(x')dx'^|(x) trFexpjzsr / A^(x")ha:"^| (0) )c • (154) 

j C/g t/ C/g 
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Here the C^’s are two infinitesimal loops centered around x ands 0 respectively, suitably defined 
on the lattice as, for example, elementary square loops. The gravitational analogue of such an 
(inhnitesimal loop) correlation was given earlier in Eqs. (1241) and (1251) . It is also understood that 
in gauge theories the inverse of the correlation length ^ corresponds to the lowest mass excitation 
in the gauge theory, the scalar glueball with mass uiq = 1/^. If the lightest scalar glueball has 
a mass of approximately m = 1750MeH (which then fixes ^ = 1/m), then in QCD is about 
eight times smaller, which gives rise to what has been described in QCD as “precocious” scaling. 
So while the two scales are quite close, they do not necessarily coincide. And the same could be 
true in gravity. 

Another important difference between gravity and QCD is that fact that in the former the 
ultraviolet cutoff still appears explicitly, hidden in the physical value of Newton’s constant G. 
There exists then a second scale ^ whose magnitude is not directly related to the value of G; 
instead it reflects how close the bare G is to the ultraviolet hxed point at Gc, and is therefore 
hne tuned in this approach (just like a mass squared term in a scalar held theory). In QCD on 
the other hand the scale ^ appears explicitly, wheras the ultraviolet cutoff is deeply hidden in the 
renormalization group relationship between the bare coupling at the cutoff scale as (A). 
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Appendix 

A Details on the parallel code 

A few details will be given here regarding the performance of the computer code used on the 
Datura Supercomputer at AEI. The latter is a high performance Infiniband cluster, primarily used 
for large scale numerical calculations in classical General Relativity. The cluster is based on Intel 
Xeon X5650 2.66 GHz processor boards which have two processors per board, each with six cores. 
Thus there are 12 cores per board, 200 nodes and 2400 cores total. The main communication and 
storage network is based on an Infiniband QDR 324-port 40Gbit/s low latency and high bandwidth 
switch. In addition, the Datura cluster is equipped with 4,800 GB of main memory, corresponding 
to 2GB per core. 

There are presently two main versions of the lattice quantum gravity code, a scalar (sequential) 
one and a parallel code. Both codes run exclusively in double precision (64 bits). The scalar code 
takes about 1021 seconds per iteration on a lattice with 32^ = 1,048,576 sites, which includes 
updates with an action containing higher derivative terms. The scalar code performance on a 
single core is measured (by counting raw floating point operations using a hardware performance 
monitor) at 3.3 GFlops. 

In the case of the parallel code, on a 16^ lattice (1,572,864 simplices) the edges emanating from 
256 sites can all be updated in parallel, with the work distributed on 256 cores. A whole lattice is 
therefore updated in 256 passes, and it takes about 0.26 seconds per iteration for a whole lattice. 

On the 32^ lattice (25,165,824 simplices) the edges emanating from 256 sites are all updated 
in parallel, with the work distributed again on 256 cores. A whole lattice is therefore updated 
here in 4096 passes, and it takes about 4.01 seconds per iteration for a whole lattice. Thus for 
this setup the parallel codes is about 255 times faster than the single core scalar code. The MPI 
communication overhead between nodes for this setup accounts for only about 1%. Using 256 cores 
the overall code performance on this machine is around 845 GFlops. 

On the 64^ lattice (402,653,184 simplices) the edges emanating from 256 sites are all updated 
in parallel, with the work distributed on 256 cores. A whole lattice is therefore updated now in 
65,536 passes. This then gives about 64 seconds per iteration. When a full complement of 1024 
cores are used instead of 256, the time for one full lattice iteration goes down to about 16 seconds. 
With a full 1024 cores used in parallel the overall code performance on this machine is around 3.4 
TFlops. If even more cores are used, the code performance should further improve provided the 
network bandwidth is commensurate. 
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